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Abstract 

We give a partonic interpretation for the deeply virtual Compton scattering (DVCS) measurements 
of the HI and ZEUS Collaborations in the small- 2:3 region in terms of generalized parton distribu- 
tions. Thereby we have a closer look at the skewness effect, parameterization of the t-dependence, 
revealing the chromomagnetic pomeron, and at a model-dependent access to the anomalous grav- 
itomagnetic moment of nucleon. We also quantify the reparameterization of generalized parton 
distributions resulting from the inclusion of radiative corrections up to next-to-next-to-leading 
order. Beyond the leading order approximation, our findings are compatible with a 'holographic' 
principle that would arise from a (broken) S0(2,l) symmetry. Utilizing our leading-order findings, 
we also perform a first model-dependent "dispersion relation" fit of HERMES and JLAB DVCS 
measurements. From that we extract the generalized parton distribution H on its cross-over line 
and predict the beam charge-spin asymmetry, measurable at COMPASS. 
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1 Introduction 



The electron/positron-proton collider experiments HI and ZEUS at the HERA ring in DESY 
improved not only the quantitative understanding of inclusive processes, e.g., by pinning down 
the small-x behavior of parton distribution functions (PDFs) jl],!^; but also led to new insights into 
the proton structure. Two decades ago, it was mostly unforeseen that at high energies the deep 
inelastic scattering (DIS) cross section steeply rises, that the proton remains intact in almost one 
third of all scattering events, and that even exclusive processes become measurable, contributing 
considerably to the total cross section. The reader may find comprehensive reviews in Refs. 

These exclusive processes, e.g., the electro- or photo-production of a vector meson or a photon, 
were extensively studied in the small-xe kinematics by the HI and ZEUS Collaborations 6N26|. 
The amplitude of the subprocesses, 

l^*\qi)p{Pi) ^ V{q2) p{P2), V = j,p, u, 0, J/^, T, (1) 

is necessarily dominated by t-channel exchanges that carry the quantum numbers of the vacuum. 
The large amount of HERA data calls for a phenomenological description and it challenges the 
theoretical understanding of the nucleon in terms of its partonic substructure. Needless to say, a 
quantitative understanding of the parton dynamics will be crucial at the frontier of exploration of 
the structure of matter at LHC 



27 



]. In this context, it is worth noting that exclusive Higgs 
production via gluon fusion is a rather clean channel 2^, 30]; however, for cross section estimates 



the gluonic content of the proton must be quantified. 

One might have hoped to master the phenomenology of such exclusive processes in the frame- 
work of Regge theory. Unfortunately, for an 'incoming' virtual photon ^-matrix theory is not 
applicable. Thus, Regge theory loses the theoretical foundation and might possibly be replaced 
by a pragmatic Regge phenomenology 3l|. Consequently, firm conclusions valid for on-shell scat- 
tering, like the one that unitarity requires that the rightmost singularity in the complex angular 
momentum plane belongs to a J = 1 exchange, might not be appropriate for off- shell processes. In 
fact, one of the lessons of Hl/ZEUS experiments is that cross sections, effectively parameterized 
as 

rise steeply, contrarily to what is implied by the pomeron (J = 1) trajectory 

ap(t) = 1 + 0.25 t/GeV^ (3) 

In addition, the effective trajectory a{t) varies with the photon virtuality = —qf- The corre- 
sponding effective Regge pole in the complex angular momentum plane might also be understood 
as a convenient implementation of cuts. 
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Inspired by QCD, a phenomenological description of cross sections at small xb has been 



achieved in terms of the color dipole model 



32 



making direct contact to the high energy 



35|. The physical picture might be set up 



approximation (BFKL) of scattering amplitudes [34, 
in the rest frame of the proton in which the highly energetic virtual photon fluctuates into a 
quark-antiquark pair. The quark-antiquark pair forms a small color dipole, spatially distributed 
in transverse direction, that interacts with the proton by a gluonic exchange. Finally, the quark- 
antiquark pair forms a meson or annihilates into a photon. The physical amplitude is given as 
convolution (with respect to longitudinal momentum fraction and transverse separation) of the 
color dipole spectral function (cross section) with the corresponding wave functions, describing 
the transition of the initial photon into a quark-antiquark pair or the quark-antiquark pair into 
the final state. 

A perturbative QCD framework, applicable for longitudinally polarized photons, is founded on 



factorization theorems 



36| . Here, in setting up the partonic space-time picture, one may prefer 



a 'brick wall' frame. The proton is viewed as a bunch of partons travelling along the ^-direction 
close to the light-cone, while the virtual photon goes in the opposite direction and may have 
zero energy. During the scattering process, e.g., a quark is knocked out by the photon, it picks 
up an antiquark and they form a meson, which travels in the direction of the photon. Another 
subprocess is when the photon knocks out a quark-antiquark pair arising from the fluctuation of 
gluonic components in the proton wave function. In this picture the quark-antiquark pair moves 
almost coUinearly and its transition into the final state is described by a collinear distribution 
amplitude. The partonic content of the proton is encoded in generalized parton distributions 
(GPDs) 37H39|. which can be interpreted as a probability amplitude to emit and absorb a parton 

ight cone or, equivalently, as a t-channel exchange 40|. For comprehensive 



that moves along the 
reviews on GPDs see 



41 



42|. 



Obviously, the two space-time pictures are related by a boost along the z-axis and in an ex- 
act QCD treatment both approaches must yield the same results. The rise of the cross section 
([2]) at large W requires a pomeron-like t-channel exchange, encoded either in the color dipole 
spectral function or in GPDs. The difference between the two approaches, e.g., for the case of 
longitudinally polarized photons in hard exclusive meson production, is mainly in the treatment 
of transverse degrees of freedom and in the view on partons, exchanged in the t-channel. This 
partonic exchange is in the color dipole model approach commonly viewed as entirely gluon dom- 
inated, while the collinear perturbative framework provides an implicit separation, defined by the 
factorization scheme, of the quark and gluonic content of the proton. The correspondence of the 
objects in the two frameworks is obvious, and the color dipole spectral function is often expressed 
in terms of so-called /c^-unintegrated GPDs or, equivalently, so-called quantum phase space dis- 
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tributions 



43| . We add that in the two-gluon exchange color dipole model one may integrate out 



the /c^-dependence in the target related part, which provides then 'hybrid' models that contain 
preasym pto tic corrections due to the transverse quark- ant iquark pair fluctuation, elaborated, e.g., 
in Refs. |4#j46|. 

In this article we restrict ourselves to DVCS at small xb, which has been measured at DESY on 
the electron/positron-proton collider HERA experiments HI and ZEUS js, ol, 7, 10|. This has been 



)reviously studied in the spirit of the aligned-jet model [47|, from high-energy/Regge perspective 



48 



50 



(LO) 



in color dipole model 



56 



58 



51I-I55I 



and in collinear factorization approach at leading order 



next-leading order (NLO) 



59 



6l|, and next-to-next- leading order (NNLO) 62|. As 



said above, the view on the partonic content of the proton varies with the approach. In DVCS 
the different points of view are obvious. In the collinear DVCS approach one starts with the 
hand-bag diagram, while in a color dipole model one would draw a t-channel gluon ladder. In 
the collinear factorization approach one usually stays with the resolution conventions which are 
set in the standard PDF analysis of deep inelastic scattering (DIS) at a given fixed order in the 
coupling. Whether one expresses now the sea quark content in terms of gluons, as in the color 
dipole approach, or stays with a large amount of sea quarks, might be considered as a matter 
of taste. However, one should bear in mind that what we shall call gluons is not the same 
object that appears in a model approaches. Let us also add that a model-dependent evaluation of 
power suppressed corrections and their interpretation as higher-twist corrections in the collinear 
factorization approach is likely an oversimplified view on the dynamics of QCD. 

We analyze the DVCS data of the HI and ZEUS Collaborations within the collinear factor- 
ization approach, going along the lines of previous work j62|. To set up non-perturbative GPD 
models, we are to some extent motivated by the Regge phenomenology and we employ its language. 
However, to avoid any misunderstanding, we enclose in quotation marks Regge-theory terms, ap- 
pearing in our modelling. For instance, "pomeron" denotes a t-channel exchange with vacuum 
quantum numbers responsible for a steep rise of cross sections, associated with an exchange of 
a colorless quark-antiquark or gluon pair. A Q^-dependence of our "pomeron trajectory" will 
be solely induced by evolution 63|. Using a least square fitting procedure, we aim for a model- 
dependent access to both quark and gluon GPDs at the LO of perturbation theory and beyond. 
There are two reasons to update our previous findings. First, the HI Collaboration provided new 
data from the HERA II run including a significantly improved measurement of the t-dependence 
of the cross section ^ and a preliminary result for the beam charge asymmetry (BCA) [10]. The 
second reason is that, in contrast to our previous ad hoc model study, we utilize here for the first 
time flexible GPD models, allowing us to describe DVCS also in the LO (hand-bag) approxima- 
tion. Such GPDs can then be used in a simple GPD description of fixed target experiments, as 
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pointed out in Ref. [64 1. 

The outline of the paper is as follows. In Sect. [2] we present the observables, relevant for our 
DVCS analysis in the small-XB kinematics. In Sect. [3] we give a short overview of popular GPD 
models that have been employed for small-xe phenomenology. Thereby, we adopt the terminology 
of quantum mechanics for both CFF and GPD models, i.e., we distinguish between GPD model 
and its representation. A model refers to a specific GPD which can be given in different equivalent 
representations, e.g., within double distribution representation, some version of conformal partial 
wave expansion, directly in some SO (3) partial wave expansion or in the light-cone wave function 
overlap representatioiij. We shall also emphasize that the claim, that at small x GPDs should 
be rigidly tied to PDFs, cannot be mathematically justified and, moreover, its partonic/physical 
content is rather speculative. We will then set up our model for integral (conformal) Mellin 
moments. Then we give a short insight into the technicalities of the (conformal) Mellin-Barnes 
integral representation of the DVCS amplitude. In Sect.lHwe present then a detailed analysis of HI 
and ZEUS DVCS data and provide a GPD interpretation. In particular, we shall discuss the failure 
of ad hoc small-x models at LO, and we shall give a detailed analysis of the skewness effect, the 
parameterization of t-dependence, the model-dependent access of the anomalous gravitomagnetic 
moment, and the reparameterization effect of radiative corrections. In Sect. [5] we include our LO 
GPD findings in a global "dispersion relation" fit to DVCS data for unpolarized proton target, 
which includes the measurements of HERMES, CLAS, and HALL A Collaborations. Thereby, we 
aim for a first model-dependent extraction of the dominant GPD H. This will allow us to predict 
the beam charge-spin asymmetry for COMPASS kinematics. Finally, we summarize and conclude. 



2 Deeply virtual Compton scattering at small xb 

In the deeply virtual electroproduction of photons both the Bethe-Heitler (BH) bremsstrahlung 
and the DVCS process contribute to the cross section. In the small-xe region the DVCS process 
dominates which allows one to extract the DVCS cross section by a subtraction procedure. There 
the integration over the azimuthal angle projects on the transverse DVCS cross section and 
guarantees that the contamination by the interference term is negligibly small. In this way, 
both the HI and ZEUS m Collaborations measured the cross section of DVCS on an 

unpolarized proton in dependence of the photon virtuality Q^, the center-of-mass energy W, and 
the momentum transfer squared t. 

^In the literature specific GPD models are often associated with representations. We have experienced that this 
indiscriminate use of language can lead to confusion. 
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On the theoretical side, the unpolarized DVCS cross section in the small xe-kinematics, 



da 



DVCS 



dt 



2c2 



Ana' 
Q4 w^ + Q? 



\n\ 



t 



4M2 



1^1 



(4) 



is primarily given in terms of two Compton form factors (CFFs), and S. Here a is the elec- 
tromagnetic fine structure constant and Mp is the proton mass. Note that we consider CFFs as 
functions of the symmetric scaling variable C,, which can also be expressed in terms of the Bjorken 
scaling variable xb: 

Both CFFs belong to the parity- and charge-even sector and, caused by an effective "pomeron 
exchange", they might rise steeply at small xb- In the coUinear factorization approach their 
dominant contribution arises from the twist-two GPDs H and E. In the cross section we 
neglected contributions of parity-odd CFFs "H and S, which are supposed to be suppressed by 



56 



0(xb), for detailqj see, e.g., Refs. 



62|. Owing to the integration over the azimuthal angle 
0, there are no interferences of twist-two with twist-three or gluon transversity contributions. 
The squares of the both latter and possible twist-four and higher contributions are all considered 
negligible. 

Furthermore, the possibility to have both electrons and positrons as probes in the HERA 
experiments allows direct access to the BH-DVCS interference term. This serves as an experi- 
mental consistency check of the aforementioned subtraction procedure and additionally allows the 
measurement of the beam charge asymmetry (BCA) 

d~^a — d~a 



A 



BCA 



(0) 



(i+cr + d a 



(5) 



reported in Ref. 10|]. This asymmetry has a more intricate azimuthal angular 0-dependence, 
particularly when one integrates over a restricted phase space. For unbinned data, it can be 
expressed in terms of the BH and DVCS amplitudes and approximated to twist-two accuracy by 
the first and third harmonics. Here and in the following we take the conventions of Ref. 56| and 
find 

BH^DVCS + ^BH^DVCS 



^bca(0) 



IT, 



BH 



+ |Ti 



DVCS 



xb- 



cos 



)+Xb 



Fi gfJe (St + 2Hr) - ^2 5Re (Ht - ^Ht) 



MT{(i>) 



(6) 



cos(3</)) , 



^Let us add that similarly to polarized parton densities the Regge phenomenology of the CFFs V. and £ might 
not be so well understood. The CFF ^5 contains also a pion pole contribution which yields a constant real part. 
However, it cannot compete with % in the DVCS cross section. Hence, for small x-b, we can safely neglect these 
CFFs in (01). 
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where Fi{t) and F2{t) are the Dirac and Pauh form factors of the proton. We could again neglect 
the "H contribution in the first harmonic. The third harmonic is induced by a photon helicity 
flip of two units. It is perturbatively tied to gluon transversity and might be contaminated by 
twist-four (quark) contribution 65|]. The normalization factors, 



y 



4M2 2-2y + y^ 



(7) 



and 
AT 



yQ 



8ynr^{2 -2y + y')Vtmin-t 



n=0 



c™ cos(ra<j 



rPi(0)P2(0)x^^c°^^ScosM) 



n=0 



moderately induce a further azimuthal angle 0-dependence. Thereby, the DVCS Fourier coeffi- 
cients cj^^*"^ depend of course on the CFFs. Although a good approximation of the normalization 
factors is obtained by restricting to the zeroth harmonics in the BH and DVCS amplitude squared 
terms. 



^DVCS 
-0 



2{2-2y + y' 



t 



4M2 



1^1 



(9) 



a residual 0-dependence is induced by the product "Pi (0)^2(0) of BH propagators. At twist-three 
level the zeroth and second harmonics of the interference term complete the azimuthal angular 
dependencies in the BCA ([6]). Note that this zeroth harmonic is dominated by the same twist-two 
CFF combination as the first harmonic, displayed in Eq. 



For the HERA kinematics we expect from the model studies in Ref. 56|], see discussion in 
Sect. 7.1 there, a moderate twist-two modulation (~ 10 — 15%) of the BCA with a small twist- 
three admixture (~ 2 — 5%), which is mainly governed by the twist-two CFF "H. We emphasize 
that the cos(30) harmonic, related to gluon transversity, is theoretically and phenomenologically 
uncharted. To reveal the dominant twist-two cos(0) harmonics of the interference term from a 
BCA measurement, one might utilize the Fourier decomposition of the charge asymmetry 



Abca{<P) =Po+Pi cos(0) + P2 cos(20) + cos(30) + 



(10) 



given as an infinite sum, and extract the dominant amplitude pi from a fit. In this way both 
twist-three and possible gluon transversity contributions diminish and we could on the theoretical 
side safely employ the approximation for the amplitude 



where denotes the zeroth harmonic of the product 'Pi{(j))V2{(f))- Obviously, the (unbinned) 
observable pi is suppressed at large and small —t and vanishes for —t — t- — tmin ~ 0. Note 
also that twist-two contributions to other harmonics might arise from the residual 0-dependence 
in the normalization factor as well. 

From cross section measurements one extracts in the first place the CFF "H, since S is kinemat- 
ically suppressed by a factor —t/AMp, see Eq. (j4]). Both CFFs enter in a different combination 
in the first harmonic of the BCA (fTTil . where ^eS is now accompanied by —F2{t)t/4Mp. A 
model study, in which S was varied with fixed H, has shown that it might be possible in a BCA 
measurement, supplemented by the DVCS cross section measurement, to access the CFF S 561 ] . 



_ To understand clearly which GPD degrees of freedom are accessible in a DVCS measurement 
.\q4\. we remind that the CFFs satisfy "dispersion relations", see, e.g., Refs. 67, 68, 66, 62 



701 1. Literally, we have to distinguish between the physical dispersion relations that follow 
general considerations and the "partonic" ones that are in one-to-one correspondence 



66 



691, 

rom 



62 



69 



70|] with the leading approximation of CFFs in the perturbative framework. For finite photon 



virtuality the physical and partonic dispersion relations differ by higher twist contributions, see, 
e.g., Refs. 6^,l70|. According to our goal, we employ the set that is consistent with the perturbative 



framework at leading twist and denote them with quotation marks. For instance, H and S have 
even signature and their subtracted "dispersion relations" might be written in the following form 



62|: 



3fte 



7r 



-PV / dx 



2x 



-3m 



x^ 



:i2) 



where we set the threshold ^th = 1 + 0(1/ Q^) equal to one. Since the real part is determined by 
the imaginary part, one mig ht choose to extract it from the experiment indirectly, by means of 
the "dispersion relation" 6J]. Then, when fitting the data, one needs to model only the imaginary 
part and the subtraction constant. It has been pointed out in Refs. 



66 



62 



3 



that the subtraction 



constant is entirely related to a GPD term that completes polynomiality. In the representation of 



Ref. 



7l| it is the so-called D-term, see also discussion in Ref. 72 1. 



The functional form for the CFFs in the small-xe region might be borrowed from the Regge 
theory. Thereby, we assume that the "Regge trajectory" is linear in t and is same for both DIS 
and DVCS. The imaginary part of CFF can then be written as: 



Q'm 



(13) 



For on-shell scattering, i.e., = 0, one would take for a{t) = a{t, Q"^ = 0) the soft pomeron 
trajectory ([3]) as the leading one. It appears in the (target) helicity conserved CFF "H and we 
shall assume that it couples also to the helicity- flip CFF £. This is related to the so-called 
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chromomagnetic pomeron, studied in Ref. 73| in the instanton approach. Interestingly, that 
study suggests that the chromomagnetic pomeron should play a dominant role in polarization 
phenomena at high-energy and so one might expect that the CFF £ is sizable. In principle, the 
pomeron might also couple to the set of transversity CFFs, which would be visible as a sizeable 
cos(30) harmonics in the BCA (|6]). 

In accordance with phenomenology and perturbative QCD the "pomeron trajectory" depends 
on the photon virtuality Q^. Note that perturbation theory states that in the double log approx- 
imation, for large and small ^, the "pomeron" in the Regge-inspired ansatz (fT3|) is replaced 
by a more intricate functional dependence 63j. Also for this reason, we understand the param- 
eterization fll3p as an effective one. Next, we have non-dominant "Reggeon exchanges", e.g., 
a{t) ~ 1/2-1- t/GeV^, which we shall neglect in the small-xe kinematics. 

Plugging the Regge-like parameterization f|T3|) of the imaginary part into the "dispersion rela- 
tion" (fT2l) . we of course recover the known Regge formula for the small-xe region: 



n 



i — cot 



ha 

6q 



if Q') + 



(14) 



Here the ellipsis stands for non-dominant contributions of "Reggeon exchanges" with typically 



a 



^'=§(t~0.2GeV2) -0.3, 



for a constancl, and for further terms that vanish as ^ — ?■ 0. In the considered kinematics the 
"pomeron trajectory" has a typical value of 

a{t ~ 0.2GeV^ ^ 4GeV2) ~ 1.1 - 1.2. 



Hence, in contrast to on-shell forward Compton scattering at large energies, where the classical 
pomeron intercept a(0) = 1 yields a vanishing real part, we encounter in DVCS with a{0) > 1 a 
positive real part that should even dominate the negative "Reggeon" contributions and subtraction 
constant 56|. This already predicts that for Hl/ZEUS kinematics the sign of the leading cos0 
harmonic in the BCA (|6]) is positive and even provides an estimate of its size, quoted above, which 



is consistent with preliminary measurements of the HI Collaboration 10|. 



■^From the viewpoint of the "dispersion relation" ([T2l) this constant, i.e., the so-cahed J = pole, contains 
contributions from the subtraction constant and further contributions induced by both "pomeron" and non-leading 
"Regge exchanges" . 
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3 Models for the flavor singlet GPDs at small x 



In the collinear factorization approach the CFFs T = {HjE} are given as convolution of a hard 
scattering part with GPDs F = {H, E}, which to LO accuracy reads 



2\ LO 



E 



q=u,d,s,- 



dx 



^ — X — iO ^ + X — iO 



(15) 



where is the fractional quark charge. Note that this is nothing else but the so-called hand-bag 
approximation, where only quarks^ are resolved in a hard Compton scattering process. Obviously, 
the imaginary part of the CFFs is given by the GPD on the cross-over line {rj = x) 



q=u,d,s,- 



el [F'^ix = e, r/ = e, t, Q') - F'^ix = -^,v = e, Q' 



(16) 



The GPD at negative momentum fraction x = — ^ corresponds to the negative antiquark con- 
tribution. The real part of the amplitude can be calculated either by means of the "dispersion 



relation" or 



Tom the convolution formula fll5l) . The equality of the two results is ensured by 



construction [68|, |62|, i.e., by the polynomiality or support properties of the GPD 66|, |62|. For 
fixed the LO perturbative approach is completely equivalent to the (approximated) "dispersion 
relation" f[T^ and the only information that can be accessed at LO is the GPD on the cross-over 
line. 

Perturbation theory also predicts the evolution of GPDs, e.g., their change on the cross-over 
line in the flavor singlet sector is 

Q2^F(x, X, t, Q') = f'^ y (1, y/x- rj = x, a,{Q))F{y. x, t, Q') , (17) 

where the flavor singlet quark (S) and gluon (G) GPDs are collected in the column vector 

[F'^{x,---)-F'^i-x,---)] , (18) 



F{x,7],t, Q' 




x,r],t, Q 



F^ix,--- 



q=u,d,s,' 



and the evolution kernel is a two-dimensional matrixo which leads to the mixing of the quark and 
gluon GPDs. F^{x, 77, ■ ■ ■ ) and F^{x, rj,- ■ ■) are antisymmetric and symmetric in x, respectively, 
and both are symmetric in 77. One realizes from the evolution equation f|T7|) that the scale change 



^ Of course, this is a simplified partonic view, since an infinite number of longitudinal gluons, radiated from 
the struck quark, are included in the GPD. The radiation of a transverse gluon is taken into account in genuine 
twist-three contributions. 

^It is common to adopt the convention in which the gluon GPD reduces in the forward limit to xG(x). In such 
a case the non-diagonal entries in the evolution matrix are accompanied by 1/?/ or 77, which are set equal to x, see 
Ref. Ol- 



io 



of the GPD on the cross-over hne is governed by its value in the outer region {x > r]). Moreover, 
the evolution kernel in the gluon-gluon channel has a "pomeron"-like singularity. Hence, as in the 
forward case, the evolution in the small-x region is driven by the gluons and can be solved in the 
double log approximation. There one finds that the "pomeron intercept" a(0) effectively incieases 
and the slope a' decreases with growing 63| . We note that the change of the residue function 
in ( |T3l) under evolution cannot be derived for the general case, since it is determined by the GPD 
in the outer region. 

Beyond LO accuracy, the imaginary part of CFFs is given by a convolution of the hard- 
scattering part with the GPDs, analogous to the evolution equation. As in DIS, the momentum 
fraction integration runs from x = ^ to x = 1. However, there is an important difference. Namely, 
the GPD itself depends on the scaling variable ^. Hence, the deconvolution cannot be performed 
and the outer GPD region can only be accessed through evolution effects, see discussion in Ref. 75 1. 

According to all this, a parameterization of the outer GPD region is sufficient for a description 
of DVCS via imaginary part of CFFs. Then, using simple Regge-motivated small-^B parameter- 
ization of imaginary part of CFFs from f[T^ . and in accordance with the LO relation flTB]) . we 
write a rather flexible GPD ansatz in the small-x region at a given input scale: 



F{x > T], T], t) = np r{^\F) (3{t\F) ; i^ = t]/x. 



(19) 



Here we factorized the GPD into the Regge part x~°'^^\ the skewness function r(r]/x\F), the 
residue function f3(t\F), and the residue rip oi the GPD in the forward limit. Hence, in this 
limit, i.e., r/ — and t — 0, both r{ri/x\F) and f3{t\F) are normalized to one. Note that 
skewness- and t-dependence might not necessarily factorize. Since present data do not allow to 
address their possible functional interplay, we stick here to the most convenient model. Also, we 
would like to add that in this GPD approach restoring of the GPD in the central region [x < rf) 
is a mathematical problem, which is solved in the short distance operator product expansion 
framework by construction 68|, |62|, |6J] . Thereby, the Mellin moments of F{x,x'd,t) are expanded 



in a Taylor series around the point i? = 0, i.e., r{'d\F) possesses certain analytic properties. 



3.1 Survey of GPD models at small x 

GPD modelling can be done in different representations. It is popular to employ either double 



distribution (DD) representation 



76 



37 



79 



40 



76l | or some version o 



sol ], adopted from the description of mesons 



conformal partial wave (PW) expansion 



84| . Mathematically, there is a one- 



to-one mapping between GPDs in two different representations. Since this mapping might be 
cumbersome, it has been sometimes only partly worked out. In these circumstances, one might 
hope to gain some new physical/partonic insight in the small-x region just by a change of the 
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representation. However, it is clear that all popular models finally lead in the small-x region to the 
GPD of the form (ITU]) , where, however, the skewness function might be rather restricted. Surely, 
the small- a: GPD properties, described in momentum fraction representation, can be equivalently 
set up in any other representation, too. 

Our effective Regge-motivated ansatz f|T9|) is supported by a diagrammatic analysis from the 
t-channel view 67|, where it was found that the "Regge trajectory" is skewness- independent; 



however, the residue function, i.e., the skewness function r('i9|F) in our ansatz ( IT9|) . depends on it. 



One arrives at the same conclusion from an s-channel view, if one adopts the ideas of Refs. |85l. |86| 
and convolutes a spectator quark model, expressed in the DD representation, with a constituent 
quark mass spectral function. Of course, in Regge-inspired DVCS models one starts with such a 
point of view, where the skewness dependence of the impact form factors has to be modelled. We 
already stated that Regge behavior is included in the popular GPD models, too. 

Such arguments and agreement with experiment, see also Sect. [H led to a broad acceptance of 
effective Regge behavior for off-shell processes, even in the absence of a firm theoretical foundation. 
Assuming this, the remaining task for GPD phenomenology is the determination of the 

i. skewness function r{i)\F) and 

ii. residue function f3{t\F). 

Certainly, one can address these questions by fitting to the experimental measurements using one's 
favorite GPD representation. The problem remains always the same, namely, to have appropriate 
functional GPD ansaetze. 

The value of the skewness function at ^9 = 1, see Eq. f|T9|) . is an important characteristic of a 
GPD model. For the GPD H, with r(^9) = r(^9|if) and r = r(l), it can be expressed as the ratio 
of the GPD at the cross-over line to the corresponding PDF, given as the GPD at 77 = for t = 0. 
Since we rely on the universality of the "pomeron trajectory", this skewness ratio is for small x 
mostly independent of x. For the quark GPD it reads 

r^Q') = ^'^"^^ s7x!q^7 ' = H''{x,v = 0,t = 0, Q') , (20) 



and for gluons we use the convention: 

^G(g2) ^ H^ix,v=^x,tj^O,Q') ^^^^^ ^ ^G(^^ = 0, t = 0, Q^) . (21) 



In an LO DVCS analysis the quark skewness fl20l) is given as the ratio of ^mn{x,t = 0, Q^)/tt 
and the transverse unpolarized DIS structure function Ft{x, Q^)/x; therefore, it can be almost 
directly measured, while the gluonic one (12T]) can be accessed only by utilizing a large lever arm 
in Q\ 
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The aligned-jet model considerations in Ref. 47| provide estimates for both the DIS structure 
function, given as imaginary part of the virtual forward Compton amplitude, and for the imaginary 
part of DVCS amplitude. In an LO approximation resulting quark skewness ratic^ turns out to 
be r^(Q^ ~ 2.5 GeV^) ~ 1. We note that this model has been used to predict the size of the 
DVCS cross section, which was afterwards experimentally confirmed, see Ref. 0]. The model was 
then generalized for the deeply virtual Com pto n amplitude with two virtual photons, providing a 
prediction for the quark skewness function |6l| : 



FMcDS 



l + ^?)l-"(0) 

2^ 



1 + 



Ml 
Ql 



In 



-12 -M2 



1 _ ^0~"'^0 ,Q 



Ml ~ 0.5 GeV^ , Ql ~ 2.5 GeV^ .(22) 



In the double distribution ODD) representation a GPD that has support in the interval —7] < 



X < 1 is expressed as follows 



37 



76|: 



F{x,r],t) 



dy 



l-z 



6{x - y - 1] z)f{y,z,t) . 



(23) 



Note that in this specific representation polynomiality is not completed for H and E. This 
artifact is not crucial for the small-x application and it can be cured in various ways, see, e.g.. 
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Refs. H 

ansatz (RDDA) 
namely. 



It is quite popular among phenomenologists to utilize Radyushkin's DD 
Here the DD at t = is factorized in a PDF and a profile function. 



h{y,z,t = 0) 



r(3/2 + 6) q{y) 



(24) 



r(i/2)r(i + 6)i-y V (i-y)'y 

The profile function is chosen to be convex and its width is controlled by the b parameter. The 
skewness function can be easily evaluated in terms of hypergeometric functions: 



r 



RDDA 



m = 2Fl 



a/2,(l + a)/2 
3/2 + 6 



' RDDA 



«-l)/2,(a)/2 



^9- 



(25) 



3/2 + 6 

Here and in the following we use a shorthand for the intercept a = a{t = 0). Setting ^ = 1 in the 
skewness functions f l25|) . the values of the skewness ratios at the cross-over line rj = x follow: 



_22.-ar(3/2 + 6)r(l + 6-a) 



'^RDDA 



r(3/2)r(2 + 26-a) ' " r(3/2)r(3 + 26 - a) ' 

Taking the originally proposed value 6 = 1, one finds for a = 1.2: 



'^RDDA ~ 1- 



'^RDDA~l-04, 



(27) 



^Thc authors define quantity R^^^ as ratio of imaginary DIS and DVCS amplitudes. In LO approximation this 
means r = 2'^(°)/R^^^ ~ 



13 



For growing b the quark skewness ratio decreases and rapidly approaches the value = 1, which 
corresponds to setting GPD equal to a t-decorated PDF. For the gluon GPD we see that the ratio 
= 1 is almost realized for 6=1. A skewness ratio r < 1 requires to give up the convexity of 
the profile function. 

It is also rather popular to expand GPDs in tree-level conformal partial waves. In some versions 
the conformal partial wave amplitudes are then additionally mapped to forward-like functions. In 
such expansions the evolution operator is diagonalized at LO. Unfortunately, such representations 
involve certain mathematical subtleties, which have been understood only in the last few years 
Originally, based on an oversimplified view on the inverse integral transformation 



91 



92 



40 



94j . it has been erroneously stated that for small x the GPD transform reduces to a PDF. Thus, 



the authors of Ref. 
the following ratio: 



94j | claim that the GPDs in the small- a; region are at t = tied to PDFs by 



_ 2"r(3/2 + a) Q _ 2^+°r(3/2 + a) 



r(3/2)r(2 + a) ' r(3/2)r(3 + a) ■ ^^^^ 

This ratio is in fact a Clebsch-Gordan coefficient, e.g., occurring in the conformal partial wave 
expansion of the product of two currents (it can be seen, e.g., by taking a "Regge pole" value 
for complex conformal spin, j = a — 1, in Eq. (1651) below). Therefore, we shall call it conformal 



ratio. Although the small-x claim, which corresponds to a specific double distribution model 95|, 
is not based on general grounds, the arguments were repeated and refined in a more recent version 
of the claim, where it was assumed that the "pomeron pole" is the rightmost singularity in the 



complex conformal spin plane 96 L Below we shall spell out, what is already (implicitly ) s aid 



in the literature, e.g., Refs. |95|, |9l|, |92|, |40|, |93|, namely, that the claim in Refs. |9J] and [96| is 
based on certain mathematical simplifications. As long as the partonic content of the conformal 
PW expansion is not understood, we consider this claim as not necessary applicable for GPD 
phenomenology. 

We would like to mention that GPDs can also be described by a t-channel SO (3) partial 
wave expansion, formulated within complex-valued angular momentum J. Of course, assuming 
effective "Regge poles" , this description is fully equivalent to our Regge-motivated ansatz f[T^ in 
momentum fraction space. The GPD representations are somehow the analogue of a Legendre and 
)ower expansion of amplitudes at high energies. We recall that both expansions are equivalent 



97( 1 and the latter became more popular in Regge phenomenology in the sixties of the last century. 
A combination of t-channel SO (3) and conformal partial wave expansion has been proposed 



79|, where the difference of conformal spin j + 2 and t-channel angular momentum J, is used as 



a discrete variable: 

p = j + l- J={0,2,4,...}, (29) 
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which is always even. This GPD representation is set up with forward-like functions Qp{z) of 
momentum fraction- like variable z and has been called "dual" parameterizatiorj^. Below we employ 
essentially the same model within a conformal partial wave expansion in terms of a Mellin-Barnes 
integral. Taking in such a model only the leading contribution, where j + 1 = J, one finds the 
conformal skewness function and ratio fl28|) . Taking into account that the singular behavior of the 
non-leading contributions Qp{z) increases by , one can obtain other values for the skewness 
ratio 931]. Hence, p cannot be a priori considered as an expansion parameter. Note that in the 

the second forward-like function 



minimal "dual" model of Ref. 



57| and in its corrected version 



Q2{z) has the same small- 2; behavior as the leading one Qq{z). Hence, this model possesses the 
conformal skewness ratio (1281) . 

The inverse transform in the "dual" parameterization has not been derived so far. Nevertheless, 
it has been exemplified within model studies, where p can be effectively employed as an e xpansion 



99 



100 |. 



parameter, that "dual" and DD representation can be effectively mapped into each other 

In the RDDA the conformal ratio fl2Sl) arises when choosing 6 = a, and, moreover, the con- 
formal skewness function coincides then with the RDDA one, given in Eq. fl2Sl) . Since the value 
a 1.2 is rather close to the originally proposed value for 6, the conformal ratios 



1.65 



1.03 



(30) 



are rather close to the RDDA ones (!27|) . In the "dual" model of Ref. [57(| both the leading Qo{z) 
and next-leading Q2{z) forward-like functions have been taken into account; however, both of 
them have the same z behavior. Hence, in this model the skewness ratio also has the conformal 
value 



and it is in the flavor singlet sector not very different from the RDDA, see also Ref. (57j 
for numerical examples. Therefore, both GPD models contradict the aligned-jet model estimate 
and so also the experiment 
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Ref. 



To convert the predictions of the aligned-jet model into the GPD language, the authors of 
6l| set the GPD at small x and low input scale p^ equal to the PDF: 



H{x, x,t = 0, pI) = q{x, pI) for small x and low pi 



(31) 



This is nothing but the RDDA model in the limit 6 — 00; practically, a large value 6 ^ 1 is 
sufficient. This ansatz implies that the skewness function is set r{ri/x) = 1 for all x. With such 
an initial condition, evolution, starting at a rather low input scale, will rapidly lead to an increase 
of the r-ratio. Thus, this GPD model fails to describe data, too. We will not go into details here. 



^Usage of term "dual" was motivated by the fact that in dual models [981 the s-channel amplitude is described 
by the t-channel exchanges. We add that this feature is more general and arises from crossing and the Sommerfeld- 
Watson transform of the ^-channel S0(3) partial wave expansion. In Regge theory /phenomenology the resummation 
of t-channel exchanges is encoded in the Regge trajectory. 
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however, we would just like to point out that the skewness function determines the evolution of 
the skewness ratio. In RDDA ansatz its functional dependence ensures that the conformal ratio 
will be approached with increasing Q^, as numerically demonstrated in Ref. 101 1. 

Let us only shortly discuss the functional form of the residual t-dependence. Historically, it is 
common in Regge phenomenology to model it by an exponential ansatz, i.e. 



f3{t) = exp {Bt} 



(32) 



In Regge framework it is clear that this is an effective description for the small —t re gion . An 



exponential i-dependence of GPDs arises in the overlap wave function representation |102| - |104 | 
within wave functions that have an exponential fall-off in the transverse momentum. It is ques- 
tionable that in s uch GPD models the GPD spectral properties, ensuring polynomiality, can be 
implemented 105| . On the other hand, power-like wave functions allow us to represent GPDs hav- 
ing both properties. With power-like wave functions one can easily evaluate DDs and they meet 
the field theory inspired view. In this framework a simple spectator quark GPD model was set 



88| and it was found that at small x the t-dependence indeed factorizes and that a power-like 



up 



behavior arises: 



t 



(33) 



Certainly, in such models Regge-like behavior is not implemented. As stated above, this can be 
achieved by convolution of such model with an appropriate constituent mass spectral function. 
However, it remains unclear to us how should the t-dependence be treated, so th at common Regge 
behavior arises and, simultaneously, general positivity constraints lod . 



Ref. 



107l | are satisfied; see 



108l | for GPD integral representations. Hence, we have no definite ansaetze for the residue 



function at smaller values of t at hand; however, a residual t-dependence as in Eq. (133!) at a (very) 
low input scale and smaller values of x looks to us rather plausible. 



3.2 Modelling of integral conformal GPD moments 

Our goal is to have a first empirical look at a more flexible GPD modelling, where we rely on the 
physical Regge-inspired picture. We shall introduce three different GPD models for the small- a; 
region, each having different skewness function. The models are defined by their conformal Mellin 
moments. This is foremostly a technical point, allowing us to employ existing stable numerical 
code that includes radiative corrections at NNLO level. Analogously to the singlet quark and 
gluon GPDs in Eq. (ITS]) , we collect their conformal moments in a two-dimensional vector 

H,{rj,t,^')= I I {v,t,fi'). (34) 
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These moments are given by the expectation values of local operators which are the lowest state in 
conformal towers; they are labelled by integral conformal spin j + 2, and they have spin projection 
j + 1. They can be defined by convolution of GPDs with Gegenbauer polynomials 

HM^imm^ir^^n'i^i' „° v^ut..,.,,'). (35) 



where polynomial indices u = {3/2,5/2} are group-theoretically determined. Since Cj{—x) = 
{—iyCj{x), our GPDs ( |T8l) have only odd j moments. For the GPD E we adopt the analogous 
decomposition. If polynomiality is completed, the integral conformal moments fl35|) are polyno- 
mials in rj of order j -|- 1. Note that the highest order terms provide the subtraction constant in 



66 



62|. 



the "dispersion relation" ([T^ 

The normalization in fl5S]) is chosen in such a way that the basic property of GPDs H, namely, 
that they reduce in the forward limit to the PDFs, translates into similarly simple integer-j 
moment-space relation 

H,{r],t,fi^)-jq^{t,fi^) = (^^^^ {t,fi^), (^^^^ {t,fi') = j^dxx' (^^^ {x,t,fi'), (36) 

where zero-skewness GPDs S and G at t = are given by the fiavor singlet quark and gluon 
PDFs, respectively. For the overall normalization we shall use the PDF momentum fraction 
averages A^^'*^, that are given by the first Mellin moments and satisfy the following sum rule 

N^in'^) + N^ijj,'^) = 1 with = Ei{t = 0, fi^) and = Gi{t = 0, fi^) . (37) 

In our modelling we are going along the lines pointed out in the momentum fraction space. 
The dominant small-x behavior arises from the "pomeron exchange" that is for zero-skewness 
GPD related to sea quark contributions defined via anti-quarks q 

S = g^^" + g™^ , q^''^ = 2q , = ^ [q-q] ■ (38) 

q=u,d,s,--- 

The difference of quark and anti-quark PDFs we denote as valence-like flavor singlet contributions. 
They are related to "Reggeon exchanges" with an intercept a ^ 1/2. Thus, they can be neglected 
for the small-XB kinematics. The standard DIS terminology fl38l) we adopt for GPDs, too. First 
we set up the zero-skewness GPD, including the t-dependence, and afterwards we skew it, i.e., 
model the skewness function. 

We start with the standard ansatz for PDF Mellin moments at a given input scale /ig 

n(^\ - Ar a:""(l-x)^ Mellintransform _ B (1 - a + J , (3 + 1) 

^^^^-^%(2-«,/3 + l) ^ ^^'"^^ i?(2-a,/3 + l) ' ^ ^ 
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for both quark and gluon PDFs. Here a > 1 is the intercept of an effective "pomeron trajectory", 
while f3 parameterizes the large x or j behavior. The normalization is given by the momentum 
fraction average A^, see Eq. (1371) . 

Then, to obtain the model for the zero-skewness GPDs fl36|) . we decorate the PDF Mellin 
moments fl39l) with t-dependence by (i) extending the "Regge intercept" a to a linear trajectory 
a{t) = a + a't, where we only introduce the leading pole, and then (ii) by multiplying with 
a residue function f3{t), having either exponential (l32l) or power-like ( l33l) functional form. In 
accordance with our t-factorized ansatz at small x, we neglect the dependence of the slope B or 
cut-off mass M on All parameters can be separately adjusted for sea quarks and gluons, even 
the effective "pomeron trajectories" might slightly differ. Thus, our ansatz for the Mellin moments 
of the zero-skewness GPDs reads 

it) - AT,,, _ ^^^^^ ^^^^ + 1) 1 - + ' ^ - ' ^ ^ 

Here the ellipsis in the quark singlet sector indicates valence-like contributions, cf. ( 138|) . The 
"pomeron pole" is written as a monopole form factor, where the cut-off mass squared rrij depends 
on j. The pole is located at 

rn^j =t =^ j = ait) - 1 • 

Note that at t = these pole contributions are included in the beta functions, appearing in 
Eqs. dinigl]). 

The skewing is achieved by relating the conformal GPD moments for given conformal spin to 
the forward Mellin moments by a linear transformation that depends on rj, e.g., by 

/ sca\ / ^sea/- N Q \ 

HM t) = r.M 1^1^ j (i) + . . . , rM = ' J"^ j . (42) 

It is required that for integral conformal spin r^^'^ir]) are even polynomials in r] of order j + 1, 
which are at 77 = normalized to one. They are moment-space analogue of the skewness functions 
r{r]/x). 

Here, we should already point out a crucial property of rjit]), given as a function of two 
variables j and 77. On the first naive glance one might expect that normalization rj(0) = 1 implies 
skewness ratio r{6 = !) = !; however, this is not the case. The inversion of Mellin transform, such 



^This is also sufficient from the pragmatical point of view because such a dependence anyway cannot be con- 



strained by fitting small-xe DVCS data [62 1 
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as (135|) . involves an integration over complex-valued j with certain weight, e.g., the GPD on its 



cross-over line is given by 40 1 



If one extends j to complex values and assumes a smooth forward limit of rj (77) and uses the pole 
value J = a — 1 in Eq. fj43l) . one is immediately led to the conclusion that the skewness ratio is 
equal to its conformal value (|28|) . e.g., for sea quarks. 



limr'^^^Cr?) = 1 r'^^i^ = D = ^"^(^/^ + _ (44) 

miw^ ^ ij r(3/2)r(2 + a) ^ ' 

This is a naive thinking, too, because assuming that the function rj(?7) for complex-j has the 
same 77 — limit as its integer-j analogue is unwarranted in the case of a singularity at 77 = 0. 
We will show in Sect. 13.31 that indeed a branch point can develop there. Precisely this feature 
ensures the equivalence of representations and will enable us to construct flexible GPD models, 
not constrained by the conformal skewness ratio ( l28l) . 

We proceed with our modelling for integer j. To find the convenient functional form for rj{ri), 
it is helpful to start by evaluating it from a known GPD at t = 0. A simple example has been 



given in Ref. 40|, where the following quark-antiquark unsymmetrized GPD has been evaluated 



from a toy DD P°^{y, z) = y °, taken to be equal to its "Regge" piece: 



1 , 1 / X 



H'°y(x,r]) = e(-7]<x)-{—^] +7] ^-7]. (45) 

1 — a 7] \1 + r] J 

At small X this toy GPD provides in the forward limit the realistic PDF behavior and its 
skewness function at small 77 is 



2^ a 



(46) 



which is reminiscent of the aligned-jet model one (l22l) . Since this expression arises from a DD 
that does not vanish on the support boundaries, it is ill-defined at = 1 for a > 1 and needs a 
regularization, e.g., by using a profile function that vanishes at the support boundaries. Note that 
for RDDA f l24p the GPD behavior at the end-point x = —77 and cross-over point x = rj is governed 
by both the Regge behavior of the PDF and the end-point behavior of the profile function. 

Nevertheless, we can use this model as a guide for modelling conformal GPD moments. The 
conformal moments of the GPD psj) can be straightforwardly calculated by means of a so-called 
beta transform of Gegenbauer polynomials into the hypergeometric function 3F2. We write the 
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resulting conformal moments in the form of a transformation ( 142|) and read off: 

(l+j-a)(2+j-a)r(3/2)r(3 + j) 



X 



2(1 -a)(2 

Tj — 1 



a] 



T{3/2 + j] 



(47) 



2r] 



-J,3 + j,2- 
2,3 - a 



a 



V 



27] 



{v -v} 



Note that for odd j these moments exist even for a = 1. The divergent behavior of the model 
( H5l) at -i? = 1 shows up in this space as convergency problem for large j. Hence, we can employ 
the skewness moments p7|) only if the zero-skewness GPD decreases faster than at j — )■ oo, 
which is satisfied for realistic values of /3 in fHOj) and fHTj) . As already said, below we shall continue 
the odd integer j conformal moments to complex- valued j. Here we only emphasize that r^°^ [rj) 
has a zero at j = a — 1 which will cancel the "pomeron pole" contribution of zero-skewness GPD. 
This obviously demonstrates that the normalization condition for odd j will not necessarily imply 
the limit (l44l) and so the conformal ratio (!28|) cann ot be taken for granted. 

In our modelling we also follow the suggestion of 
in SO (3) t-channel partial waves 62 1 
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79| and expand conformal GPD moments 



(48) 



labelled by the angular momentum J, where rj = —1/ cos 6 is expressed in terms of the t-channel 
center-of-mass scattering anglqj. The PW amplitude -Fjj+i with angular momentum J = j + 1 
is nothing but the Mellin moment of the zero-skewness GPD, while the PW amplitudes with J 
smaller than j + 1 are suppressed by a factor r]^ (p = j + 1 — J). The SO (3) partial waves djr, 
normalized to one at r/ = 0, are the crossed version of Wigner's reduced rotation matrices (ig a-a' 
4l[ |. For the t-channel helicity conserved 'electric' (A = A' = 1/2) and helicity flip 'magnetic' 
(A = -A' = 1/2) GPD moments 

t 



4M2 ' 



and Hi + E 



respectively, they are expressed by Legendre (Gegenbauer with index u = 1/2) polynomials, 

t 



4(r/) oc r]^dlo{l/v) 

_ r(i/2)r(j+i) 

2-^r(l/2 + J) 



for T = n 



4M2 



£ 



(49) 



(50) 



^We are here in the first place interested in the assignment of quantum numbers. Thus, we neglect some 
corrections which die out either in the Bjorken limit or by setting the proton velocity in the c.o.m. frame to one. 
These corrections appear in the partial waves and their amplitudes. However, after completing the Sommerfeld- 
Watson transform they will partly cancel each other in the small- a; region. Hence, it is for our application more 
convenient to neglect them from the beginning. 
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with Jjnin = and Gegenbauer polynomials with index v = 3/2, 

for T = n + S, (51) 

r(i/2)r(j) ,,0/2 



2-^r(j+ 1/2) 

with Jmin = 2. Although we here consider integral J, it is already appropriate to think in terms of 
the Regge language, providing us a guideline for the resummation of fHHj) . The quantum numbers 
associated with our conformal GPD moments are associated with "pomeron" and parity- and 



charge conjugation-even "Reggeon exchanges" llOl . l4l| . where the leading ones belong to the 
/-meson trajectory, generically given as a(t) ~ 1/2 + t/GeV^. Each term in the sum (148|1 should 
therefore have a leading pole at J = a(t) rather than at j = a(t) — 1. Note that the former 
scenario, well established from the physical point of view, yields a non-trivial skewness ratio, 
while the latter, more formal and vague scenario, yields the conformal ratio (144 p . Polynomiality is 
completed, i.e., the "magnetic" and "electric" GPD moments are polynomials of order j + 1 and 
J — 1, respectively, where j is odd. Terms proportional to ri^~^^ contribute to the so-called D-term, 
which finally builds up the subtraction constant in the "dispersion relation" . To this constant not 
only the J = a- meson contributes, but rather all exchanges with J < j + 1. Still, in our small-xe 
DVCS application the "pomeron trajectory" dominates these "Reggeon exchanges" , including the 
subtraction constant, and so we shall hereafter ignore them in our small-x GPD modelling. 

As one realizes from Eq. ( HHjl . for t = the proper partial waves for Hj are Legendre poly- 
nomials (I50p . while Ej can be represented as a sum of —Hj and an addendum, expanded with 
respect to Gegenbauer polynomials f l5T]) . For the sake of simplicity we take these assignments also 
for t ^ 0. Moreover, in the small-?] kinematics the different SO (3) partial waves have (for complex 
valued j) the same asymptotic behavior. Hence, in the following we will relax the requirement of 
the complete polynomiality. 

The SO (3) PW expansion (HHI) can be written in the form of the transformation (H2|) . One 
immediately reads off that the SO (3) PW expansion of the skewness matrix is given by 

rMF)= E -^V'^'-'d^iv)- (52) 

even 

We now construct three different models for the skewness dependence, which will be employed for 
fitting of small-XB DVCS data. Showing also the corresponding "dual" model we have: 

i. a leading SO (3) partial wave (1-PW), i.e., minimalist "dual" [p = 0) model, 

ii. a leading and next-leading S0(3) partial wave (nl-PW), i.e., the minimal "dual" (p = 0,2) 
model, and 
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in. a model-dependent resummation of SO (3) partial waves (S-PW). 



i. l-PW model: Taking in the expansion fl48p only the leading SO (3) PW J = j + 1 into account 

1 /2 

and representing the Legendre polynomials by hypergeometric functions, we write this 

model in the form: 



1 



7] — 1 



27] 



1 
1 



(53) 



where the forward moments Qj are specified in Eqs. (HUj) and f l4ip . This is the most restrictive 
model, and its skewness ratio is fixed by the conformal value as claimed in Refs. [9J, |96 |. 



ii. nl-PW model: To build a flexible model we include also the next-leading partial wave amplitude 
with J = j — 1. For simplicity, we express this PW in terms of the leading one, where, however, 
the conformal spin is shifted by two units: 

s^)' ^^^^ 

The ?7^-proportional term, having a pole at J = j — 1 = a(t), is crucial for adjustment of the 
normalization of the DVCS amplitude, controlled by the entries s^'^^ and s*^ in the matrix S. For 
complex conformal spin the pole is at j = 2 + a — 1 within the constraint 3fJej > 3 and might thus 
be denoted as spurious. 

Hi. 'L-PW model: It turns out that the sum (HHj) of SO (3) partial waves with physically-motivated 
effective "Regge poles", 

J - a 

yields after analytic continuation in j the same small-?] asymptotic as implemented in our toy 
skewness moments P7|) . see below Eq. f l69|) . To have a flexible skewness ratio, we take the l-PW 
model and add the difference of the toy model P7|) and the l-PW onj^: 

iff (r/,t) = Hf{rj,t) + S,irj,t\a) q^{t) , (55) 

where the entries of the skewness matrix 

^ / X I s^'^^,(r/|asea(t)) \ 

' V s^SMo^oit)) 



^"Note that this second term wih die out in the forward hmit, and so it is invisible in the zero-skewness GPD. 
Of course, such a trick can be used in any representation. 
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are defined in terms of the skewness moments (147|) and (153|1 : 

For convenience, we included an additional factor r(7+a) /r(7+a(t)) that makes the t-dependence 
of skewness effect more fiat. 



A few comments are in order. Skewness parameters are denoted by the same symbols, s^'^'^ 
and s^, in both nl-PW and S-PW models; however, their normalizations are not related in an 
obvious way. Still, for both models positive (negative) values of these parameters imply increase 
(decrease) of the normalization of the DVCS amplitude at t = 0, compared to the conformal 
ratio situation (1251) of the 1-PW model. The existence of both these flexible models explicitly 

n 

shows that the claim of Ref. [9J] about conformal ratio cannot be derived from a conformal PW 
expansion (or equivalently from the integral transformation 78|). The assumption of Ref. 96l ]. 
namely, that singularities in the complex conformal j plane with Kej > a — 1 (where a ~ 1) are 
absent would exclude the chosen representation, including the constraint 3fJej > 3, for the nl-PW 
model. To weaken the constraint, we will change the Mellin-Barnes integral representation for 
CFFs (or GPDs) by a variable shift j — ?■ j' — 2, as it is done below in Eq. f l67|) . Then the spurious 
pole will be moved to the "right" position J = j' + 1 = a{t). This is in accordance with our 
assumption that the leading poles are associated with the angular momentum J rather than with 
the conformal spin j. Moreover, to have a mathematically consistent representation of the S-PW 
model we will perform below the analytic continuation by means of Carlson's theorem and then 
the spurious right half-plane poles in the complex j plane are absent anyway in the conformal 
moments Utilizing the small-?] expansion, we will also illuminate that spurious poles should 
be considered as an artifact of this expansion. 

Finally, we specify the models for the GPD E. Since in the small-// region the SO (3) partial 
waves fl50|) and fl5T]) behave rather smoothly and they can be safely approximated by one, we do 
not need to be careful about specific choices of expansion polynomials. Hence, for simplicity, we 
can assume that the conformal moments of Ej are proportional to those of Hj: 

where the ellipsis stands for valence contributions. The normalization of the GPD moment Ej is 



for j = 1 given by the anomalous gravitomagnetic moments, see, e.g., Refs. lllNll3| . of flavor 
singlet quarks and gluons 
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Ji's decompositioiJ^ of the proton spin [115| reads then 



(59) 



As a consequence of both momentum (157|) and angular momentum conservation ( I59p . the sum 
rule 



B^ + B*- 







(60) 



states that the anomalous gravitomagnetic moment of the nucleon vanishes. Note, however, that 
the residue of Ej at the "pomeron pole" j = a — 1 depends besides the normalization B also on 
the functional form with respect to j. 

The flavor decomposition of the anomalous gravitomagnetic moment is an important phe- 
nomenological goal, which we will address below in Sect. (14. 5 p within the ansatz ( 158|) in an ad 
hoc model-dependent manner. Let us set up a few scenarios. Generic thoughts within simplified 
models, utilizing large x counting rules, isospin symmetry, and the values of the nucleon mag- 
netic moment, suggest that the sum of valence u and d contributions provide an almost vanishing 
valence- like anomalous gravitomagnetic moment 116|. With this taken literally, we would have 
the scenario: 



B' 



On the other hand, lattice simulations 1171 Ill8| suggest the scenario 



(61) 



(62) 



which would suggest that in our model the chromomagnetic "pomeron" should be absent. However, 



disconnected contributions, which are relatec 
account in present lattice measurements 117 



;o a gluonic t-channel exchange, are not taken into 
118|, and one might be tempted to reinterpret the 



lattice results as valence dominated, i.e., B = 0, supporting our generic thoughts (16T]) . 



3.3 Numerical evaluation of CFFs 

As motivated above, we approximate the CFFs J-' = {H,S} in the small-xe region by their flavor 
singlet part 



- E 

Nf ^ 



(63) 



q=u,d,s- 



^^A overview on the common spin decomposition schemes can be found in Sect. 5 of Ref. 114| . 



24 



where the (averaged) fractional squared charge Cg is for Nj active quarks. In the factorization 
formula for the CFFs we will set both renormalization and factorization scales equal to the photon 
virtuality Q. Instead of the momentum fraction representation, we utilize the conformal partial 
wave expansion for complex-valued conformal spin. The flavor singlet part of our signature even 
CFFs is then evaluated via a Mellin-Barnes integral 

J 1 -L ran I 

2 



1 

2i 



i + tan 



(64) 



The coefficient functions C and the evolution operator E are evaluated in perturbation theory. 
For instance, in LO (hand-bag) approximation the hard scattering amplitude reads 



. LO 2^-+^r(j + 5/2) 

r(3/2)r(j + 3) ^ ' ' 



(65) 



while the evolution operator is the same as in unpolarized DIS. The radiative corrections to the 



62| for the standard 



coefficient functions C and the evolution operator E are presente in Ref . 
minimal subtraction (MS) scheme in NLO and for a special conformal scheme (CS) up to NNLO. 
We add that in the momentum fraction representation the corresponding expressions are known 
for the MS scheme in NLO approximation^^!. 

A few comments are in order. In the CS scheme conformal symmetry is manifest, up to the 
breaking by the trace anomaly, also beyond LO. Thus, the NNLO corrections can be obtained 
by conformal mapping from the corresponding forward quantities, elaborated in DIS; for detai ls 
we refer to Ref. 



62| and references therein (for integral conformal spin see also Refs. 



123 



m). 



In the CS scheme the convolution [C ® E]^. reduces to a simple multiplication, while in the MS 
scheme an additional summation over the conformal spin must be included. 

In the evolution operator we resum only the leading logs and perturbatively expand the non- 
leading ones. In the small-^ region the perturbative expansion of the universal evolution operator 
is getting unstable, so that the difference to an evolution operator in which also non-leading 
logs are resummed is of some numerical importance. This difference is implicitly absorbed by a 
reparameterization of sea quark and gluon PDFs or GPDs at the given input scale. 

The continuation of coefficient functions and anomalous dimensions from integral to complex- 
valued j is, as in DIS, straightforward. This is also the case for zero-skewness GPDs, already given 



^^For clarity, the perturbative quantities which were denoted in 62| as C and £, are renamed in the present 
paper into C and E. 

^^Since even at NLO a diagrammati c calculat ion of evolution kernels has not been performed yet, except for the 
two simplest ones in the quark sector 119l4l22l |. we do not expect that the NNLO results in the MS scheme will 
become available in near future. 
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in Eqs. (140| HTj) as ratios of F functions. The continuation of the skewness moments, as polynomials 
of 1] is more intricate. It is required that these moments obey certain bounds for j — )■ oo with 
|a'rg(j)| < 7r/2, see Ref. 401]. Then Carlson's theorem tells us that their analytic continuation is 
unique. The conformal moments ( 153|) of our 1-PW model, expressed by hypergeometric functions 
2-Fi, and their descendants, i.e., the hypergeometric functions 3F2 appearing in the skewness 
moments (1471) with argument (77 — l)/2?7 < 0, satisfy the required condition. However, their rj- 
symmetrized counterparts, with argument {1 + ri)/2r] > 0, are problematic. Such functions are 
now evaluated on the branch cut in the complex rj plane. The correct treatment is to take the 
discontinuity over the cut, and to add a term that restores polynomiality for odd j: 



3-^2 



3-^2 



-j, 3 + j, 2 - a 
2,3 - « 



V 



1+2 

2r, 



It 



27] 
3^2 



(66) 



1+2 

2r? 



2i sin(7rj) 



2t] 



2-a 



T{l+j)T{3-a)T{l + a + j] 

r(3 + j)r(«)r(3-a + j) 



In the continuation procedure for odd j moments the factor {—riY is replaced by —rjK The Mellin 
moments, obtained by this continuation procedure, do not possess singularities for JJej > a — 1 
with a ~ 1 and are unproblematic for a < 2. It has been already numerically checked for a < 1 in 



Ref. j40| that our toy model in the conformal Mellin space is equivalent to its momentum fraction 
space representation. We convinced ourselves that for t] ^ ^ this holds true also for a ~ 1.2. As 
already mentioned, the singular behavior of the toy GPD at the cross-over line is absent in our 
S-PW model. 

For the numerical evaluation of the Mellin-Barnes integral (164|) we choose as integration path 
a straight line segment, parallel to the imaginary axis. Intercept c of this line and the real axis is 
taken to be c ~ 0.35, so that the contour lies to the right of the "pomeron pole" and to the left 
of the poles of tangent function in f l6^ . For the nl-PW model the polynomiality condition of the 
?7^-proportional term in Eq. f l5^ requires us to take a separate Mellin-Barnes integral where the 
integration path is shifted by two units to the right, i.e., the intercept of the integration path is 
now c + 2. Shifting now the integration variable, j ^ j — 2, we combine both pieces into a single 
Mellin-Barnes integral, schematically written as 



1 

2i 



c+ico 



tan I — 



+2 ' 



■Hf 



(67) 



for our nl-PW model. Obviously, the ?7^-suppressed term in integral conformal moments con- 
tributes to the leading Regge behavior and influences the value of the residue function. We 
emphasize again that the spurious pole in the nl-PW model is now moved to the "right" position 
J = j + 1 = a and that the new integration variable can now be viewed as the (shifted) t-channel 
angular momentum, i.e., j = J — 1 . 
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For the small-^ kinematics, we are interested in, we can speed up numerics by expanding 
the conformal GPD moments in the vicinity of 77 = 0. For 1-PW model fl53p . using well-known 



expansion of Legendre functions 125l | one gets 




^\V) = ( 1 [1 + 0{v')] + 0{v'^+') H%,t) = q^{t) . (68) 



With our choice for the integration path, i.e., 3?ej = c > 0, the 0{rf^^'^) term is negligibly small 
and the conformal GPD moments reduce to the Mellin moments of zero-skewness GPDs. The 
same approximation we can use in Mellin-Barnes integral fl^Tj) for the nl-PW model. The small-?] 
expansion of the skewness moments fH7|) for the S-PW model, for complex valued j, obtained 
within the prescription (!66|) . was worked out in Ref. 



40|: 



(toy)_^ ^r/y+i-r(l/2)r(l+j) r(l-a)r(l + a+j) 2[l + 0(r/2)] ^ 



' \2) r(3/2 + j) 4"r(a)r(l-a + j) l + tan(f )tan(^) ^' 



(69) 



This expression has a series of "conformal sibling poles" in the left half-plane at j = —a — 1 — k 
with k = 0, 1, 2, ■ ■ ■ . The series of spurious poles in the right half-plane at j = a + 1 + 2k with 
k = 0, 1, 2, ■ ■ ■ are absent in the exact result, and they appear here as an artifact of the small-?] 
expansion. The function rj(?/) has also a zero at j = a — 1 that cancels the "pomeron pole" 
in the Mellin moment; however, the corresponding CFF still grows as at small C,- We have 
numerically checked that the exact and approximate result s ag ree well even for larger values of 



7] ~ 0.2; an example for a = 1/2 was presented in Ref. [40|]. In contrast to "Regge poles", 
the "conformal siblings" will move to the right with increasing —t. Fortunately, as long as the 
condition 

-*<^ (70) 

a 

is satisfied, it is guaranteed that they all lie to the left of the lowest integral value j = 1. All of 
the 'conformal' trajectories contribute to the residuum of the leading "Regge pole" at J = a{t) 
and this provides the possibility to adjust the normalization of CFFs. 

Plugging the approximations fl68|) and fl69l) into the S'-function fl57j) of the S-PW model, 
provides us with the expression convenient for numerical treatment 

s.ivHt)) = - (Hy^''"^'^ v^r(i+j) r(i-a(t))r(i + a(t)+j) 2r(7 + a)/r(7 + «(t)) ^ ^^^^ 



27 r(3/2 + j) 4-Wr(a(t))r(l-a(t)+j) i + tan(f)tan(zif) 

Note that the factor ^J+i"" for t] = ^ is partially cancelled in the Mellin-Barnes integral ( 167|) by 
which leads to the behavior of CFFs. Hence, the Regge behavior in conformal space is 
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not entirely related to a j = a — 1 pole. This is natural from our point of view that the "Regge" 
behavior is associated with a pole in the complex J plane, and, consequently, the residue is then 
expanded in terms of conformal PWs. 

Finally, we have numerically checked that the small-?] approximation of the skewness moments 
works rather well, i.e., for ^ < 0.01 ~ 0.1) we find a per mil (below 1%) effect for the S-PW 
model. The resulting CFFs satisfy also the "dispersion relation" ( IT2|) for < ^ < 0.01 ~ 0.1) 
on the per mil (4%) level. As expected, the medium- and large-^ regions are unimportant in the 
dispersion integral for the real part at small ^. The small-?] approximation of the SO (3) PWs has 
been often considered, e.g., in Ref. 62|, and is even more unproblematic than the approximation 
of the E-PW model. 



4 GPD interpretation of DVCS data in collider kinematics 

We present now a detailed model study of the small-xe DVCS measurements of the HI and ZEUS 
Collaborations 

mm 

by means of least square fits at LO, NLO in the MS and CS, and at 
NNLO in the CS scheme. The fitting procedure is described in Sect. 14.11 Thereby, we aim for a 
partonic or GPD interpretation; however, the nature of the DVCS process precludes a coverage of 
the whole space of possible GPD models and, therefore, we are unable to determine the subspace 
of models that correctly describe the DVCS data. Instead, we use three representative models 
defined in Sect. 13.21 in an attempt to illuminate all facets of the problem. In Sect. 14.21 we give 
specific emphasis to LO fits, since the outcome might be conveniently employed in studies of fixed 
target experiments. GPD reparameterization effects due to radiative corrections are considered in 
Sect. 14.31 In Sect. 14.41 we extract the transverse distribution of partons from our fits. Finally, we 
address in Sect. 14.51 the question of whether measurements of beam charge asymmetry allow one 
to access the GPD E and, consequently, to get a handle on the chromomagnetic "pomeron" and 
the anomalous gravitomagnetic moment. The lessons from the fits are listed in Sect. 14.61 



4.1 Fitting strategies and parameters 

Our numerical studies are based on the twist-two approximation of the DVCS cross section 



and BCA ([6]) without gluon transversity, the models from Sect. 13. 2[ the Mellin-Barnes integral 
fl64p . and on the least square fitting routine MINUIT 126|. We fit to the DVCS cross section 
and to the DIS F2 structure function data, where we add statistical and systematical errors in 
quadrature and look for the values of GPD model parameters that minimize the function. For 
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instance, for the DVCS cross section a we have 

^cxp _ ^the^2 



= E rA^,^2 ^ ' = \/(^sys^d)2 + (5stat^d)2 , (72) 

where d enumerates the experimental data points. The quahty of fits is assessed by the value 
divided by the number of degrees of freedom (d.o.f.) and by probability of such and larger values. 
As bad fits we consider those where x^ /d.o.i is considerably larger than one, or, more precisely, 
where its probability is smaller than 0.2. We will also separately present the x^ values contributed 
by t-, W- and Q^-dependence. Since there is no one-to-one correspondence between individual fit 
parameters and the particular kinematic parameters of the CFF, we cannot determine the exact 
number of the d.o.f. Therefore, we simply give then x^ P^r number of data points (n.o.p.). The 
difference is small for large number of points; so this should provide good idea about the quality 
of the description of dependence on the individual variables. An error analysis of the resulting 
parameters is an intricate task, which we postpone. 

We set the input scale, at which the GPD models and parameters are specified, to 

= 4 GeV^ . (73) 

We consider then the charm quark as massless and set Nf = 4 and = 5/18. Numerical effects 
of b quarks, known to be small in DIS, will be neglected her e. A lternatively, the production of 



heavy quarks might be described within perturbation theory 127j |. worked out at NLO j9l|. The 



flavor scheme choices we consider to be foremostly related to different partonic interpretations. 
The value of the QCD running coupling is specified as 



as{Qo = V23GeV) 

271 



( 0.0606 ' 




f ^° 1 


0.0518 


> for < 


NLO 


[ 0.0488 ^ 




NNLO J 



(74) 



and corresponds to its phenomenological value ajMzo = 91.18 GeV) = 0.114 at the reference scale 



Mzo, using the standard evolution prescription 128| . However, we perform forward evolution over 



the intervals of interest keeping a fixed number of four active quarks. 

Let us recall that our GPD H models are parameterized by the averaged momentum fraction 
A^, "pomeron trajectory" a{t), large-j parameter /3, skewness parameter s, and one parameter 
that controls the residue function /3(t), e.g., the S-PW model reads 



for I = {sea, G}. In our previous studies 62| we relied on the ad hoc 1-PW model (s^ = 0) and 



performed a simultaneous fit to DIS and DVCS data, where the momentum sum rule was not 
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imposed. We utilize now this constraint to fix the gluon normahzation, 

= 1- N'"'' - A^™' , A^™i = 0.4 , (76) 

where the momentum fraction for valence quarks A^™' is compatible with values obtained from 
standard PDF fits. Our fiexible GPD models could be fitted just to DVCS data, where PDF 
normalization and "Regge intercepts" 

N""", a""", (77) 

at the input scale would be taken from some of the standard parameterization^- However, 
instead of this, we mostly perform two-step fits. First we fit to the DIS data to extract the PDF 
parameters fl77|) . There the (3 parameters for sea quark and gluon PDFs are taken as 

^sea _ = 6, (78) 

which is slightly larger than their canonical values from counting rules, accounting thus for the 
increase with resolution scale. Inclusion of DIS fit in our procedure is convenient because it ensures 
identical scheme conventions and approximations for the evolution operator in DIS and DVCS. 
The DIS data used come from HI data set 129| for the DIS structure function F2 in the region 



2.5 GeV^ < g2 < 90 GeV^ . 

The choice of these 85 data points is sufficient for our studies and at the input scale (I73p our 
resulting PDFs are consistent with standard parameterizations, as demonstrated at LO and NLO 
accuracy in Fig.[Tl Note that the PDFs in the regions 10~^ ^ x 10~^ and 10~^ < x (mainly) arise 
from extrapolation within our functional pomeron-like ansatz and that we use our conventions. 



e.g., perturbative NLO corrections are consequently expanded in powers of [62 1. 

The remaining free model parameters are the skewness parameters, the cut-off masses and 
power behavior in the residue function ( !33|) . [or exponential slope parameters in ( !32|) ]. and the 
"Regge slopes" 

s^"% s^; M^"% p^"^ [E^"^], M^, [B^]; a'^ , (79) 

which control the normalization, t-dependence, and the shrinkage of the diffractive forward peak 
in the DVCS amplitude, respectively. These parameters are to be determined by fits to DVCS 
data in the second step of the procedure. However, it turns out that from our fits we cannot 
simultaneously pin down the cut-off masses and powers p in the ansaetze fl83l) . so we fix the latter. 
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Due to the small-XB collider kinematics, we neglect the "Reggeon" contributions to the flavor singlet part. 
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Figure 1: Parton densities at LO (up) and NLO (down) are shown for quarks (left) and gluons 
(right) at the input scale Qq = 4GeV^ (solid). For comparison we also show Alekhin's 



para- 



meterization 130| with errors (vertically dashed band) and some other standard PDFs 13lL Il32 | 
(dotted, dash-dotted). 

Large —t counting rules would suggest to take for sea quarks p^^^ = 4 and for gluons = 3. To 
have the possibility of direct comparison of the fitting results with the characteristic size of the 
proton, given by the cut-off mass in the dipole parameterization of the Sachs form factors, we 
choose a dipole ansatz: 

^sea =pG^2. (80) 

Moreover, since gluons are not directly accessible and a full separation of quarks and gluons 
through the evolution effects is not yet possible, we must further reduce the parameter space. 
In our previous ad hoc model studies we learned that fitting routine tends to use the gluonic t- 
slope to adjust the normalization of the total DVCS amplitude. Here we get rid of this unwanted 
complication by fixing the gluonic cut-off mass (or slope) and taking its value from the analysis 
of elastic J/ip production. According to Ref. 4^, this process is dominated by the two-gluon t- 
channel exchange and the charm quark mass already provides an internal hard scale that translates 
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to an effective J/ip photoproduction scale of Q^^j ~ 3GeV^. Therefore, one may even view the 



J /ip photoproduction as a hard process with access to the gluon GPD 27|. From the differential 
cross section measurements 23l-l26| the "pomeron trajectory" (here defined with respect to W 
rather than ^) and the residual t-slope were extracted by a fit with an exponential ansatz 26|: 

a{t) = 1.224 ± 0.010 ± 0.012 + (0.164 ± 0.028 ± 0.030) t/GeV^ , 

b = 2B= (4.630 ± 0.060^°:?i) /GeV^ . (81) 

Also, the measurements strongly favor an exponential t-dependence of the differential cross sec- 
tion, where the t-slope b is quoted in Eq. ([H]) for W = 90 GeV. Still, at the GPD leveQ we 
consider it more natural to take a power-like form for the residual t-dependence, as argued in 
Sect. 13. 1[ Therefore, the dipole ansatz will be used for the majority of our fits, with some fits 
using exponential ansatz performed for comparison. The value of the "Regge slope" parameter a' 
is not fully pinned down. The measurements of the ZEUS Collaboration 24| gi ve a slightly lower 
mean value as quoted in Eq. f l8ip and the electroproduction measurements 25|, l26|] indicate that 
a' might be compatible with zero at a larger resolution scale. The photoproduction fit (ISTj) serves 
us to fix the slope of the gluonic residue function, while for the "Regge slope" at the input scale 
fl73|l we will consider two values. To maximize model differences, we choose the combinations, 

{M^ = VoTzGeV, a'c = O.lS/GeV^} and {B^ = 2.32/GeV^ a'^ = 0} , (82) 

for dipole and exponential ansatz, respectively. Here, the values of M*^ and B'^ correspond to 
each other for o'q = 0. 



Finally, we equate the quark and gluon "Regge slope" parameters, a'^^ 



a 



G- 



parameter sets for the DVCS fits read 

{s'''', s^, M^""} , with fixed 



/(L7GeV, a^ea = "g = 0.15/GeV^ 



Then our 



^3) 



for the dipole, and 



s^, B"""} , with fixed 5^ = 2.32/GeV^ a^. 



for exponential ansatz. 



The small- xr DVCS cross section measurements of the HI and ZEUS Collaborations are pub- 
lished in Refs. js, 9| and 0], respectively. The kinematics covers the intervals 

3 GeV^ < < 80 GeV^ 45 GeV < IV < 145 GeV , 0.1 GeV^ < -t < 0.8 GeV^ 



^^Whether the measurement of the longitudinal part of the J/ip electroproduction cross section, known to much 
smaller accuracy t han the photoproduction one [26[, is then consistent with the coUinear factorization approach, 
elaborated in Ref. 



133| to NLO, is to our best knowledge not yet investigated. 



32 




Figure 2: Two-step fit at LO to DVCS [8|,S,3 and DIS [l2g data with the S-PW model and the dipole 
ansatz (|83p: (a) differential DVCS cross section versus —t for three values of Q'^ at W = 82 GeV 0], (b) 
total DVCS cross section versus Q'^ at W = 82 GeV (circles, solid) @, ^ and at W = 89 GeV (triangles, 
dash-dotted) 0], (c) total DVCS cross section versus W at the same values as on (a) [g |, an d (d) DIS 
structure function F2{xb, Q^) versus for xr = {8 • 10-^,3.2 • 10~^ 1.3 • 10^2,5 • 10-"^} |l2a |. 



We do not include in our fits one ZEUS data selo, HI data versus xb, and mo del- dependent 
extractions of skewness ratio and t-slope given by those Collaborations. Note that the more 



recent HI measurement j^, providing a larger data set on the t-dependence, was not used in our 
previous GPD study [g^. Altogether, the HI and ZEUS measurements provide us with 101 DVCS 
data points. 

We give a LO fitting example for the S-PW model in Fig. [2l Here the first panel shows the 
t-dependence, the second and third the W- and Q^-dependence, respectively, and the fourth one 
our fit to the DIS structure function F2. 



^^The excluded data is for the cross section, integrated over t, versus W for fixed = 9.6 GeV^. Including this 
set and the analogous one from HI would unavoidably increase the value of our fit. These two data sets are 
mutually compatible; however, we give here preference to HI data set, since it possesses a smoother H^-dependence. 
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4.2 Leading order fits to DVCS cross section measurements 



To simplify our analysis of DVCS cross section measurements we will neglect the helicity non- 
conserved CFF £. According to the formula ([1]) for the DVCS cross section, this is justified if the 
condition 



-t 



holds true. Since the average value of —t is ~ 0.2 GeV^, there will be at the most a few percent 
contamination by £ as long as the modulus of this CFF is comparable to \H\. In the case that 
\£\ turns out to be considerably larger than l^l, our findings, or better to say our interpretation, 
will become questionable. 



In Sect. 13. II we already mentioned that a reasonable description of DVCS data at 
in the small-XB region could not be achieved in previous ad hoc GPD model studies 
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O accuracy 



62 



58|. On 



the other hand, one expects that this should be possible using more flexible GPD models. These 



model 


as 


xVd.o.f. DIS 


xVd.o.f. DVCS 


X?/n.o.p. 


X^/n.o.p. 


X^z/n.o.p. 


1, dipole 


LO 


49.7/82 


280./100 


181. /56 


63.6/29 


36.2/16 


1, exp. 


LO 


49.7/82 


316./100 


192. /56 


79./29 


44.9/16 


nl, dipole 


LO 


49.7/82 


95.9/98 


53.2/56 


27./29 


15.8/16 


nl, exp. 


LO 


49.7/82 


97.9/98 


49.1/56 


31.2/29 


17.7/16 


S, dipole 


LO 


49.7/82 


101. /98 


57.7/56 


27.4/29 


16./16 


S, exp. 


LO 


49.7/82 


102./98 


51. /56 


32.3/29 


18.6/16 


1, dipole 


LO 


321. /182 


189./56 


51.1/29 


27.9/16 



Table 1: values and their individual contributions, coming from experimental data on t-, W-, and 
Q^-dependence, for various models with a residual dipole (j83|) and exponential ([8^ t-dependence for 
two-step fits (first DIS, then DVCS; first six rows), and for simultaneous fit (DIS-I-DVCS; last row). 
Boldface numbers indicate bad fits, as defined in Sect. 14.11 

statements are quantified in Tabled], where for various models we list the total DIS and DVCS 
over the number of d.o.f., as well as the partial contributions to from various subsets of data, 
over the corresponding n.o.p. Our ad hoc model (1-PW) is with x^/d.o.f. ^ 3 highly disfavored at 
LO accuracy in both fitting strategies (two-step and simultaneous), while flexible (nl- and S-PW) 
GPD models correctly describe both DIS and DVCS data (for them we display only results of 
two-step fits). For both flexible models, and with either dipole (1551) or exponential flM|) residual 
t-dependence, we have x^/d.o.f. ~ 1, and partial values indicate a good description for the t-, 
and Q^-dependence. The values of our fitting parameters are listed in Table [2l 
In Section l4.2.1l we reveal the reason for the failure of ad hoc models. Then, in Section l4.2.2l we 
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model 




yysca 


a [(}) 


{ 1{ /rsca\2 
\M ) 




a 


„G 
S 


^sca 


r^cff 













[GeV J 








[Gov J 


[GeV J 




1, dipole 


LO 


0.152 


1.158 


0.062 




1.247 




33. 


5.7 


0.19 


1, exp. 


LO 


0.152 


1.158 






1.247 




29. 


5.1 


0.23 


Ill, dipole 


LO 


0.152 


1.158 


0.48 


-0.15 


1.247 


-0.81 


4.8 


5.5 


0.13 


nl, exp. 


LO 


0.152 


1.158 




-0.18 


1.247 


-0.86 


3.1 


5.8 


0.14 


S, dipole 


LO 


0.152 


1.158 


0.42 


-11. 


1.247 


-32. 


5.4 


5.5 


0.14 


S, exp. 


LO 


0.152 


1.158 




-13. 


1.247 


-34. 


3.1 


5.8 


0.15 



Table 2: Model parameters, as obtained by two-step fits from Table [H together with quark GPD H 
t-slope at X = 10-3 and = 4GeV2, OFF % t-slope (USD h"^ ai W = 82 GeV and = lOGeV^, 
both in GeV"^, and resulting BCA For the dipole and exponential (exp.) t-dependence the fixed 
variables are given in Eq. (|83p and (j84p . respectively. Again, boldface numbers arise from bad fits. 



provide some insight into skewness- and t-dependence and their cross-talk for our flexible models. 
4.2.1 The failure of the small-x conformal skewness ratio at LO 

A whole class of GPD models has the same small-x behavior as the 1-PW model, characterized 
by their skewness ratio being equal to its conformal value f l28|) . As explained in Sect. 13. H for 
the "pomeron" case, of interest here, we can also include the "dual" model of Ref. 58| and the 
original RDDA (6 = 1 ^ a) in this class. We explained in Sect. 13.21 that, contrarily to the claim 
of Refs. 9^, |96| , the conformal ratio fl25]) cannot be a general GPD property. It is important to 
clarify the phenomenological status of this small-x claim. We shall now have a closer look at its 
failure by utilizing the 1-PW model. 

The 1-PW model implies a normalization of the CFF that generally overshoots the experimental 
data, which is also manifested in the large skewness effect f l28p . i.e., r ~ 1.6. To compensate for 
the too large normalization in the fitting process the t-slope gets increased via a very low cut-off 
mass parameter [M^'^'^Y ~ 0.05 GeV^ or a very large slope parameter B^^'"^ ^ 30/GeV^, see Table 
[2J As a consequence, the t-dependence of the cross section, with x^/n.o.p. ~ 3, is particularly 



badly described. A simultaneous DIS/DVCS fit also does not help, see Ref. 62| or last row of 
Table [H We add that in such a fitting strategy the structure function F2 is well described, where 
the gluon PDF comes out slightly softer than in a two-step fit. 

We exemplify now that the failure of the 1-PW model cannot be cured by modifying the t- 
dependence. From the second line of Tabled] we see that a purely factorized t-dependent ansatz, 
namely, the exponential one fl8^ . is with x^/d.o.f. ~ 3.2 even more disfavored than the dipole 
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ansatz (183!) . Such an exponential ansatz was also employed in the "dual" model of Ref. 57|. Our 
poor description of the DVCS data contradicts the statement of Ref. 57| that such a model is 
compatible with the DVCS data of the HERA I run . We repeated the analysis of Ref. |57| 
taking a pure Regge ansatz with no residual t-dependence, but with large "Regge slope" instead: 



O.g/GeV^ and 



Or 



O.S/GeV^ . 



^5) 



and our statement that 1-PW 



96| ) are disfavored at LO [62 1 



One might expect that such a choice lowers the normalization of the CFF, 

T^yi « r(3/2 + .W)r(2 + .) 

W^ii.t) r(3/2 + a)r(2 + a(t)) ' '^^ ^' 

in such a way that DVCS data are better described. In contrast to this expectation, we have found 
total disagreement with the data, e.g., x^/d-o.f. = 2100/101. Fortunately, the inconsistency of 
our findings with those of Ref. 57| has been resolved in Ref. [58 
models (or minimalist "dual" model and the small-x claim [OJ, 
holds true for a pure Regge ansatz f lS^ . too. Note that the value of a'^ as large as in Eq. (15^ is 
already excluded by the HI and ZEUS electroproduction measurements of vector mesons, which 
is dominated by two-gluon exchanges. 

We also recall that in a previous LO investigation 6l|, performed in momentum fraction space, 
both the quark and gluonic r-ratios were taken to be one at the input scale. What happens then 
is that the large amount of gluons, with their constant skewness function, driv es th e quark GPD 
to rapidly approach its conformal skewness ratio with increasing Q^, see Ref. 10 1| . The failure 
of this model shows that the skewness ratio for gluons is also not equal to the conformal value 

^ 1. To investigate this some more, we used the 1-PW model for gluons, where the skewness 
ratio is close to one, see Eq. (!28|) . and allowed a flexible parameterization of the sea quark GPD, 
e.g., the E-PW model with dipole ansatz. We indeed obtained a bad fit 

xVd.o.f. = 278/99 , x?/n.o.p. = 168/56 , x^/n.o.p. = 70/29 , xh/n.o.^. = 40/16 , (86) 

verifying our expectations. 

It is popular in LO model descriptions of hard exclusive vector meson production 



45 



46l.ll34. 



1351 ] to reduce the amount of gluons by taking the NLO PDF parametrization rather the LO one. 



since NLO gluon PDFs are by a factor of two or so smalle r than LO ones, see error bands in 
Fig. m^b). This recipe has been also used in Refs. 



134 



135j to conclude that the exclusive J/\l/ 
production is phenomenologically consistent with the LO small-x claim of Ref. [oj] : 



R 



G 



^7) 



In our opinion this finding does not necessarily support the small-x claim, obtained from a LO 
analysis. The recipe is only justified in the phenomenological context of modelling the hard 
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exclusive production of vector mesons 45l. l46l . Il34l . Il35| : however, in a simultaneous description of 
DIS and DVCS data this recipe is obviously in conflict with the counting scheme of the collinear 
factorization approach. Moreover, if we take at the input scale NLO PDFs instead of LO ones in 
the model of the previous paragraph, we found that although 

xVd.o.f. = 154/99 , x'/n.o.p. = 95/56 , Xw/^-O-V- = 38/29 , Xg^/n.o.p. = 21/16 (88) 

improve, compared to Eq. fl88l) . the DVCS fits are still unacceptable. 



4.2.2 Skewness ratio, £-dependence, and their cross-talk 

Our flexible models resolve the normalization problem and fits are unproblematic in both two-step 
and simultaneously fitting strategies, see Tabled] and Fig. [2J As expected, we find negative values 
of the s-parameters, listed in Table [2l This can be also viewed as the reduction of the large 
skewness effect of the 1-PW model. We recall that the definition of these s-parameters does not 
directly allow to read off from them the size of the skewness effect, or to compare their size in two 
models. 

We quantify the skewness effect in the momentum fraction space by the r-ratio fl20|) . An 
analogous quantity, emphasizing the physical rather than partonic aspect, can be defined as the 
ratio of the imaginary parts of amplitudes for DVCS and forward Compton scattering, which can 
be expressed in terms of the differential DVCS cross section at t = and the DIS cross section 
<^t{i*P X) for transversally polarized photon exchange, respectively. Assuming an exponential 
t-dependence, the skewness effect is revealed by utilizing the total DVCS cross section (9|: 



R{W, Q2 



v/167raDvcs&(Q')/(l + P') lo H'''''{x, r, = x,t = 0,Q' 



Here the t-slope is extracted from a fit with a' 



3, 



X=2x/(l+x) 



19) 



b{Q'^) = A[l - B\og{Qy2GeV^)] , A = 6.98 ± 0.54 GeV-^ 5 = 0.12 ± 0.03 , (90) 

and the ratio of real to imaginary part of the DVCS amplitude might be set to p = — cot(a(Q^)7r/2). 
At LO this i?-ratio can also be expressed in terms of sea quark GPD, where the momentum fraction 
for the PDF (or zero-skewness GPD at t = 0) is X = 2x/(l + x). The relation between skewness 
ratios fl2Ql) and fl89l) . considered now as a function of x rather than W, follows from their definitions: 



R(x, Q2) = 2"(2'V(Q2) for small 



X , 



(91) 



where we set q^'^^{x, Q'^)/q^'^^{X, Q^) = 2°*^^^\ In the kinematical region considered here a ~ 1.2 
and so i? ~ 2r. 
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Figure 3: (a) skewness ratio R (j89p for W = 82 GeV^ compared to HI data and (b) quark skewness 
ratio r (j20p for x = W~^. Model parameters, obtained from fits as in Fig. [2l are specified in the main 
text and in Table [2] for 1-PW (dotted), nl-PW (dot-daslied) , S-PW (solid) with dipole ansatz and S-PW 
with exponential ^-dependence (dashed). 

In panels (a) and (b) of Fig. [3] we display R and r for fixed W and for fixed x, respectively, 



as functions of Q^. Needless to say, the 1-PW model, obeying the small- a; claim 9J, |96| of 
conformal ratio for quarks, is in conflict with experimental measurements (dotted curves). The 
experimentally measured i?-ratio is well reproduced by both the nl-PW (dot-dashed) and the S- 
PW (solid) models within a dipole ansatz (jHS]), which are barely distinguishable. One sees that 
the i?-ratio is slightly larger than 2, while the r-ratio is generically r ^ 1, with their ratio as 
expected from Eq. (19T]) . This experimental value of i? ~ 2 is sometimes interpreted as a large 
skewness effect; however, from the GPD perspective, i.e., r ~ 1, it is much more appropriate to 
consider it as a manifestation of a small or zero skewness effect. 

We also display the result for a S-PW model with exponential f l84p t-dependence (dashed). 
Compared to the dipole ansatz fl83l) . it provides only a slightly smaller skewness ratio, i.e., the 
modulus of its negative s-parameters is slightly larger. In other words, the normalization of CFF 

at t = is larger for CFF with a dipole t-dependence in order to compensate for initial faster 
decrease with t. This can be also observed by comparing the parameters B^^^ and s in the third 
and fourth (or fifth and sixth) rows of Table [2l Thus, the quite drastic correlation of skewness 
effect and t-dependence for the disfavored 1-PW model, observed in the previous section, appears 
also for flexible GPD models, but in a much milder form. Note that the CFF 

(t)^ 0.2 GeV^Q2) 
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Figure 4: Gluonic skewness ratio r*^ (j2ip versus for x = 10 (a) and skewness function r{9) (b). 
Models are the same as in Fig. [3l 

at the mean value {t) ~ 0.2 GeV^, and the effective slope over the accessible t-interval, see Eq. f l92|) 
below, are described equally good with both ansaetze. 

Let us also have a closer look at the gluonic skewness effect. The large negative skewness 
parameters s*^ in Table |2] indicate that the gluonic r-ratio is much smaller than one. This is seen 
in Fig. m^a) for both flexible models, which are again barely distinguishable, and whose gluon 
skewness is even negative at the input scale. If the skewness ratio at the input scale is constrained 
to be positive, we could not obtain acceptable LO fits. Although the GPD on the cross-over line 
f] = X has no probabilistic interpretation, and negative values are thus not forbidden, a gluonic 
GPD model with a zero is suspicious. We consider this zero as an artifact related to an improper 
modelling of the skewness function r{-&), defined in Eq. f|T9|) . at the initial scale. The fact that 
with increasing the gluonic r-ratios of both models approach in the same way the conformal 
ratio vq ^ 1, reveals that their skewness functions are almost the same, see Fig. IH^b), although 
their conformal moments look quite different. We add that the same turns out to be true for the 
quark GPD models. Hence, finally one realizes that the S-PW model is effectively equivalent to 
the nl-PW one. 

Within our models, the t-dependence is well described at LO accuracy, too, and we cannot 
discriminate between the two types of functional dependence. As argued in Sect. 13. ![ we prefer 
the dipole ansatz (183|) with a small value of a'^^,^ = = 0.15. In panel (a) of Fig. |2] a fitting 
example is shown and the values, listed in Table [H confirm a good description of the DVCS 
data set. 

To convince the reader in a more obvious way that our dipole GPD parameterization ( 183|) 
is compatible with the exponential fit ( 190|) of the HI Collaboration, we evaluate the effective 
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Figure 5: Effective t-slope ([92]) versus for W = ^2 GeV (a) and versus W for = ig GeV^ (b), 
compared to HI data 0]. Models are the same as in Fig. [3l 

exponential slope 

Thereby, we consider the experimentally accessible interval, i.e., we set 

tl = -0.1 GeV^ and = -0.8 GeV^ 

and evaluate b'^^ from our GPD models that are fitted only to the cross section measurements. 

In Fig. |5] we show the effective slope b'^^ versus for W = 82 GeV and versus W for = 
10 GeV^. One clearly realizes that the solid and dash-dot-dotted curves, which result from the 
S- and nl-PW model with a dipole ansatz f l83|) . respectively, describe the experimental HI data 
very well. In the left panel the evolution of b'^^ follows the experimental data and so the 
perturbative evolution of these GPD models is fully compatible with the measurements. Also the 
flatness of the H^-dependence of the data, indicating the absence of a shrinkage of the diffractive 
forward peak, is well reproduced with the input values a'^^^ = = 0.15 GeV^ , see Fig. |5]^b) 
Consequently, one cannot conclude that the DVCS measurements indicate that a' is zero at a 
lower scale. 

The quark dipole cut-off mass parameter M^'^'^ ^ 0.67 GeV in both of our models is rather 
similar and somewhat smaller than the fixed gluonic cut-off mass ^ 0.85 GeV. They essentially 
coincide with the characteristic scale of the nucleon, as it appears in the dipole parameterization 
of its electromagnetic form factors. Note that the dipole masses and the resulting quark slope 
B^^^ at the input scale obviously differ 10 — 15% between the two flexible GPD models, compare 
third and fifth rows in Table [2l This essentially reflects the residual t-dependence of the skewness 
effect at the input scale and establishes a cross-talk between t-dependence and skewness. 
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Alternatively, we can also fit the experimental data set by an exponential t-dependence ([HI 
with a'^^^ = = 0. As already mentioned, the quality of fits are equally good, compare third 
(fifth) with forth (sixth) line in Tabled] Our exponential t-parametrization resembles the one used 
by the HI Collaboration ^ in the fit fl90|) of the differential DVCS cross section measurements. 
Indeed, at the input scale Q = 2 GeV our LO S-PW fit gives b = 6.2/GeV^ which is consistent 
with the HI value ([90]), b = (6.4 ± 0.5)/GeVl As in the HI fit ([90]), the evolution to higher 
scales decreases the slope parameter. However, in our GPD parameterization this reduction is 
slightly weaker than in the HI fit, displayed below in Fig. [101 together with the NLO prediction 
of the E-PW. Nevertheless, we conclude that the decoration of the exponential t-slope parameter 
in GPDs by an additional dependence as in Refs. 60|, |57[, |58| is a redundant complication, 
which spoils the perturbative evolution. The exponential ansatz f l32p describes the measurements 
well, too, where the small offset between GPD description and HI measurements is a bit larger 
than for the dipole ansatz, see also HI fit below in Fig. [10] 

Let us finally have a closer look at the "Regge slope" parameters. A small value of < 
0.25/GeV^ can be deduced from J/\& data. However, one might wonder whether a large value of 
'^sea ~ 1/GeV^, as considered in Ref. 101| and utilized in Ref. 57|, plays some role for Hl/ZEUS 



kinematics. From the Regge theory point of view, see discussion in Sect. 13. 2[ such a large slope 
parameter belongs to "Reggeon exchanges" , that are suppressed relatively to the "pomeron" one 
by C,^ with A > 1/2. Therefore, we expect that our fits should disfavor a large a' value for the 
quarks. 

To investigate this, we have performed a series of fits with the S-PW model, accompanied by 
the dipole ansatz, with different fixed value of a'^^^ parameter. Probability of the of these fits 
is plotted in Fig [H] One notices that the preferred value is indeed around a'^^^ ~ 0.15 GeV~^, 
while values beyond 0.5 GeV~^ are excluded at 95% confidence level. This confirms our model 
assumption that "Reggeon" contributions are not dominant in the Hl/ZEUS kinematics. Note 
also that the choice Oq = 0, shown as dashed curve, is disfavored in the dipole ansatz. 

4.3 Beyond leading order fits 

We shall now pursue the model dependence and GPD reparameterization in DVCS fits beyond 
LO. At NLO the gluons become a part of the hard-scattering process, and they may induce large 
radiative corrections. The inclusion of radiative corrections in a DVCS fit is compensated by the 
reparameterization of the GPDs at the initial scale. Within a flexible GPD parameterization we 
have good fits beyond LO in both the MS and CS schemes. In Table [3] we list the values for 
the dipole ansatz. 



It has been demonstrated in Ref. 



62| that beyond LO a good description of HERA I DVCS 
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Figure 6: Probability of (such and larger) 'x^ obtained by fitting the S-PW model and the dipole ansatz 
with Oq = 0.15/GeV^ (solid) and = (dashed), with respect to the ' 'Regge slope" parameter aggai 
which is held fixed dm'ing fitting procedure. 



model 




xVd.o.f DIS 


xVd.o.f DVCS 


X?/n.o.p 


X^/n.o.p 


Xga/n.O.p 


1 


NLO(MS) 


71.6/82 


148./100 


77.6/56 


36.8/29 


33.9/16 


1 


NLO(CS) 


71.6/82 


105./100 


62.9/56 


25.1/29 


17./16 


nl 


NLO(MS) 


71.6/82 


102./98 


60.2/56 


23.9/29 


17.5/16 


nl 


NLO(CS) 


71.6/82 


104./98 


61.4/56 


24.9/29 


18.1/16 


S 


NLO(MS) 


71.6/82 


101. /98 


60./56 


23.9/29 


17.5/16 


S 


NLO(CS) 


71.6/82 


104./98 


61.5/56 


24.9/29 


18.1/16 



Table 3: values as in Table[T]of NLO DIS and DVCS fits for various models in the MS and CS scheme. 

data [s, ~\ was possible also within the 1-PW model, where, however, the dipole cut-off masses were 
used to adjust the normalization of the DVCS cross section, see Sect. 7.2 and Fig. 16 of Ref. 62 1. 
After including the HERA II DVCS data with many more data points for the t-dependence, 
it is, however, not necessarily true that the 1-PW model still works. In particular for the MS 
scheme is significantly large, see first line in Table El while in the CS scheme such a model fit 
is acceptable within the dipole ansatz. In contrast, within an exponential ansatz the fits in the 
MS scheme are acceptable, however, in the CS they are now disfavored: 

r 50/56 
1 108/56 

where in both cases Xvv'/^-O-P- Xgz/n.o.p. values are fine. 

These findings illustrate the intricate interplay of the functional form of a given ad hoc ansatz 
and radiative corrections. If one would have relied on the claim of Ref. 9J] about the conformal 



xVd.o.f. 



r 95/100 
\ 155/100 



Xt/n.o.p. 



for 

CS 
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skewness ratio, one might have wondered whether the outcome of the NLO fitting favors the 1-PW 
model with the exponential ansatz in the MS scheme or with the dipole ansatz in the CS scheme, 
see, e.g., the discussion in Ref. 57|. We should recall here that the difference between the MS and 
CS schemes is entirely related to 'non-diagonal' rotation effects^ that vanish in the 77 = case. 
In other words, this scheme rotation is nothing but a skewness effect and we observed in fact an 
interplay of skewness and t-dependence, discussed in the previous section in context of flexible 
models. So we emphasize again that an 'overinterpretation' of fit results within an ad hoc GPD 
model is inappropriate. 
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Figure 7: Quark skewness r (120p (a) and gkion skewness re (|2ip (b) versus for 2; = 10 ^, extracted 
from a two-step fit. Dash-dot-dotted (green) line is 1-PW model at NLO(MS), while other lines are for 
S-PW model at LO (dotted, black), NLO(MS) (dash-dotted, green), NLO(CS) (dashed, blue) and NNLO 
(solid, red). In all models we employ the dipole ansatz 



The discussion of the previous paragraph becomes superfluous, if one utilizes flexible GPD 
models. Then we can provide also in the MS scheme good two-step fits and we can reveal the 
skewness effect. The extracted NLO fitting parameters for the MS and CS schemes are listed in 
Table m for the dipole ansatz. If one goes from LO to NLO, the most drastic changes appear in 
the skewness parameters. They mutate from large negative values to moderate positive ones in 
the MS scheme. This qualitative jump is also illustrated in Fig. [TJ compare LO fit (dotted curve) 
with the others, where we show the skewness ratios for x = 10'^ versus Q^. (This ratio for a fixed 

is rather flat over the experimentally accessible interval x G {lO^'^, 10"^}.) One realizes in the 
left panel that the S-PW quark model at NLO MS (dot-dashed) overshoots now the conformal 



^^Obviously, this rotation is numerically significant in the DVCS kinematics, i.e., 77 = ^. This shows that the 
statement of Ref. namely, that such 77-proportional effects arc negligibly small and so the LO claim [9^ about 
the conformal skewness ratio remains valid at NLO, is amiss, too. The differences between both schemes have also 
been studied in Ref. [62 1. 
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model 


Us 


i V 




[IVI ) 


S 


a [u) 


S 


Dsca 
D 





T3r^ A 


1 


NLO(MS) 


0.168 


1.128 


0.71 




1.099 




3.5 


5.0 


0.10 


1 


NLO(CS) 


0.168 


1.128 


0.57 




1.099 




4.2 


5.7 


0.09 


111 


"\TT /^/tv /ro\ 

NLO(MS) 


0.168 


1.128 


0.59 


0.04 


1.099 


0.02 


4.0 


5.6 


0.09 


111 


NLO(CS) 


0.168 


1.128 


0.58 


-0.01 


1.099 


-0.01 


4.1 


5.6 


0.09 


S 


NLO(MS) 


0.168 


1.128 


0.60 


3.10 


1.099 


1.10 


4.0 


5.7 


0.09 


S 


NLO(CS) 


0.168 


1.128 


0.58 


-0.42 


1.099 


-0.58 


4.1 


5.6 


0.09 



Table 4: Model parameters as specified in Table [2] obtained by fits from Table [3l 



model as 


xVd.o.fDIS xVd.o.fDVCS x?/n.o.p xl^/^-O-V x's^/n.o.p 


1 NNLO(CS) 
111 NNLO(CS) 
S NNLO(CS) 


81. /82 80.1/100 49.2/56 19.1/29 11.9/16 
81. /82 78.8/98 46.7/56 18.9/29 13.2/16 
81. /82 78.8/98 46.7/56 18.9/29 13.2/16 



Table 5: values as in Table [T] for NNLO DIS and DVCS fits within various models in the CS scheme. 



ratio fl28l) . i.e., r^^'^ ^ 1.6. In the right panel it is illustrated that the gluon ratio at NLO matches 
the conformal ratio fl28|l . i.e., ^ 1. We recall that the gluon PDF, extracted from DIS, is in 
LO approximation twice as large as in NLO, see Fig. (H^b), and so one might consider r'^ ~ 1/2 
as a realistic LO skewness ratio. 

In the CS scheme the skewness parameters s^^^'^ are small and negative for the dipole ansatz 
(183|) and so all models are compatible with the conformal ratio (!28|) . obtained at LO. This is 
displayed for the S-PW model by the dashed curves in Fig. [71 The same holds true if we include 
in this scheme NNLO corrections, shown as solid curve. Table states that all models provide 
similar and good values. The values of the s-parameters, listed in Table El are small and 
negative and with a slightly larger modulus than at NLO for the CS scheme, cf . Table IH Hence at 



model as 


jysca c^seaj^Q^ {M'"'^)^ s"^^ a^(0) 5*^^^ b""^ BCA 


1-PW NNLO(CS) 
nl-PW NNLO(CS) 
S-PW NNLO(CS) 


0.172 1.125 0.56 1.104 4.2 5.6 0.12 
0.172 1.125 0.57 -0.01 1.104 -0.04 4.2 5.6 0.12 
0.172 1.125 0.57 -0.89 1.104 -1.80 4.2 5.6 0.12 



Table 6: Model parameters as specified in Table [2] obtained by fits from Table [5l 



NNLO in the CS scheme, same as at NLO, non-leading S0(3) partial waves do not play essential 
role. The differences between the three models are so small that we obtain the same values for 
the effective and partonic t-slope parameters. 

The reparameterization of sea quark and gluon GPDs can be directly read off from Fig. [HI 
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Figure 8: Sea quark (a) and gluon (b) GPDs within dipole ansatz ( l33i) as they result from a LO 
(dotted), NLO MS scheme (dash-dotted) and CS scheme (dashed) as well as NNLO CS scheme 
(solid) fit. 



Here the GPDs on the cross-over line are shown for the S-PW model with dipole ansatz as they 
arise from our fits at various orders and schemes. It is again obvious that the most drastic effect 
appears if we go from the LO (dotted) to NLO description. Quark GPDs at NLO, compared to 
the LO ones, are enhanced by a factor of two or so while the LO gluons suffers from our model 
artifact. The sizable reparameterization effects for quarks of about 100% or even more are in view 



of the quoted corrections in Ref. |62l) I surprising. Naively, we would have expected as in DIS 



moderate correction on the level of 20% or so. The very large NLO corrections in the quark sector 
might be connected to the artifacts of our gluon model, discussed above, and we would here not 
exclude the possibility that a fully flexible GPD model possess in the quark sector only moderate 
NLO corrections. Consequently, for such a model the quark skewness ratio must then be ~ 1 
also beyond LO, i.e., much smaller then the conformal ratio r^'^^ ^ 1.6. The difference at NLO 
between the MS (dash-dotted) and CS (dashed) scheme are clearly visible for the sea quark GPD, 
while this skewness-induced effect is tiny for gluons. This simply reflects the properties of the 
off-diagonal scheme transformation, which is set by conventions. As already observed in Ref. 62|] 
within the 1-PW model and a simultaneous flt, the NNLO corrections lead only to a slight change 
of the parameters, obtained at NLO in the CS scheme, compare the dashed and solid curves in 
Figs. [8] and also [7] or the corresponding entries in Table H] with Table [61 

One might wonder whether the observation that the CS (N)NLO skewness ratio approaches the 
LO conformal value fl2S]) is a deflnite model independent feature or simply an accident. Naively, 
one might expect that a reduction of gluonic skewness ratio will also decrease the quark one, since 
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There a given 1-PW model has been employed to estimate the radiative corrections for the CFF H. 
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they contribute with a different sign to the CFF. In this way one might reahze a quark skewness 
ratio which would be closer to the LO findings. However, taking a more negative value of s*^, like 
that occurring in LO fits, provides bad NLO fits. As in our unsuccessful attempt in the previous 
section (to have a positive gluon GPD at the cross-over line at LO), we consider also this finding 
as a model artifact. Namely, our models do not allow adjustment of the scale dependence of the 
skewness ratio in a flexible way. 



4.4 Transverse distribution of partons 

We would now like to deliver a partonic interpretation of our analysis of the t-dependence of the 
DVCS cross section. We stick to the probabilistic interpretation of zero-skewness GPDs in the 
infinite momentum frame 
of 



1361 . In this frame the proton might be viewed as a disc with a radius 



4-^lnFi(t) ^0.6fm, (93) 
at t=o 

arising from the parameterization of the Dirac form factor in terms of the common dipole ex- 
pressions for the Sachs form factors within a cut-off mass of 0.71 GeV^ We add that the electric 



charge radius \/Jr^ = 0.875 fm of the proton, quoted in the Review of Particle Physics 137 |. is 
about ten percent larger than that from the dipole parameterization. Hence, also our proton disc 
radius (!93|) would then increase to about 0.66 fm. 



The transverse width of parton distribution, i.e., the average distance y (6^) of the struck 
parton from the proton center, is directly given by the t-slope of the zero-skewness GPD 

(x, Q2) = 4^ In H{x, = 0, t, Q') . (94) 
at t=o 

Let us first consider the LO interpretation, where we take the S-PW model as it is specified in 
Table [2J As discussed in Sect. 14. 2. 2^ the skewness effect will talk back to the t-slope parameters 
on the 10% — 15% level, which implies a < 7% model uncertainty for the transverse width. The 
rigidity, still present in our flexible models, does not allow us to reveal from LO fits to the DVCS 
data the transverse distribution of gluons. Therefore, we take the parameters as before for the 
dipole 0831) and exponential 0841) t-dependence, where the gluon t-slope parameters are fixed from 
the J/ip photoproduction. 

The value of the transverse width, obtained from fitted GPDs, also depends on the functional 
form of the residual t-dependence, which defines the extrapolation from t he a ccessible t interval 
to t = 0. For the exponential ansatz our value for gluons is quoted in Ref. 138l | and coincide with 
the proton disc radius of 0.6 fm. For a' = 0, the dipole ansatz would provide a 10% larger value 
than the exponential one. Also a non-vanishing value of a' leads to an increase of the transverse 
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Figure 9: Transverse width Vlr of sea quark (a) and gluon (b) distributions versus x for = 4GeV^ 
(thin) and = lOGeV^ (thick), resulting from the LO fit of S-PW model with dipole (solid) and 
exponential (dashed) t-dependence ansatz, specified in Eqs. ([83|) and ([8^ . respectively. 

width. Taking a' = 0.15 in exponential ansatz, the gluonic transverse width increases, e.g., to 
^ 0.72 fm for x = 10^"^, coming rather close to the one in our dipole ansatz (1331) with ^0.77 fm. 
We might conclude that even for an exponential t-dependence, the experimental uncertainties of 
the a' measurement in J/ tp production do not exclude a gluon transverse width at small x that is 
larger than the proton disc radius fl93l) . 

The gluon transverse width, obtained at the input scale, is depicted in the right panel of Fig. M 
as thin solid (dipole) and thin dashed (exponential) curves. Our LO fit states then that the quark 
transverse width at the input scale is slightly larger than for gluons, where the characteristic 
x-dependence, related to the partonic shrinkage effect, appears. We infer that our LO quark 
interpretation strongly reflects our assumptions and that we cannot discriminate, e.g., for x ~ 10"^ 
and = 4GeV^, between 



0.9 ] „ f 0-8 1 , , r dipole ([83 



(62) ^ <^ ■ Km, \ m ■ Um ioT i ^-L^ ^ ansatz. (95) 

\ /sea 1 0.7/ ' V \ /g \o.6J \ exponential dHU) J ^ ^ 

That at the input scale (thin curves) gluons are more centralized than sea quarks has been also 
found in Ref. 63(| by a flne-tuning procedure. Thereby, the latter were dynamically generated 
within the double log approximation at LO from gluons with a (soft) pomeron trajectory ([3]). 
We emphasize that there the differences between a tripole and an exponential ansatz were hardly 
visible in impact space for distances smaller than 1.5 fm and that the larger value of the transverse 
width for a tripole ansatz arose due to its long tail in impact space, see illustrative examples in 
Ref. [631. 

At larger resolution scale, e.g., = lOGeV^, the partonic shrinkage effect is practically 
washed out, see the flatness of the thick sohd curves in Fig. [91 This indicates that a' rapidly 
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approaches zero with growing Q^. Also the initial differences between the transverse widths of 
quarks and gluons diminish by a slight decrease of the former and increase of the latter. However, 
the 'typical' value associated with the particular t-dependence ansatz is robust. Finally, our model 
findings for = lo GeV^ 



0.9 ) „ „ r dipole (jHS]) 



9) ^ \ m K { \im ioi { ^ > ansatz . (96) 

V\ 'g I 0.65 J \ exponential ' ^ ' 



are compatible with those of Ref. [63[, where, e.g., a value of ~ 0.85 fm was quoted for a tripole 
ansatz. 

We would also like to emphasize that for the exponential ansatz the gluon transverse width is 
at the input scale the same as quoted in Ref. 27|. It perturbatively evolves to a slightly larger 
value with increasing scale. In Fig. |9] we also extrapolate the transverse widths obtained from 
fitted GPDs to larger x G [10~^, 10"^]. This extrapolation is not supported by the interpretation 
of J/ip data from fixed target experiments in Ref. j^^], where within a dipole ansatz a value of 
^ 0.53 fm for x ~ 0.1 was quoted. Interestingly, t he in crease of this value by about 20% — 30% to 
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We add that in realistic GPD models the 



~ 0.6 fm has been explained with chiral dynamics 
t-dependence dies out at x — 1, e.g., seen in lattice simulations of Ref. 118|, and so the partons 
are entirely concentrated in the center of the proton. This feature can be simply implemented in 
our models by decorating the cut-off mass (or the exponential slope parameter) with j-dependence. 

Next we shall show that our LO findings with fixed gluon slope parameters are only moderately 
influenced by perturbative corrections and scheme conventions, or even by the release of the gluon 
slope parameters at NLO (or NNLO). Note that this was not the case in our previous investigations, 
where we employed the 1-PW model (62|. We mainly work in the MS scheme to NLO accuracy, 
which is appropriate for a possible global GPD analysis of deeply exclusive electroproduction 
processes in momentum fraction space, and we use the S-PW model. We might have employed 
the nl-PW model as well, however, beyond LO there is only very little model uncertainty induced 
by the skewness dependence, see B^^^ values in Tables E] and O 

We first demonstrate in Fig. [10] that at NLO the description of the measured t-slope is for 
fixed gluon slope parameters fully analogous to our LO findings in Fig. |5l Again, the effective 
slope parameter (p2|) is evaluated from the outcome of our fits to the DVCS cross section with the 
dipole (dot-dashed) and exponential (dotted) ansatz, where the NLO (MS) parameters are given 
in Table HI Comparing the corresponding curves in both figures one can barely see a difference. 
We also display the HI fit fl90|) as solid curve. Obviously, only if experimental errors could be 
very drastically reduced, one might be able to discriminate between the three curves — for a more 
detailed discussion see Sect. 14.2.21 

A simultaneous release of all four i-related parameters (i^'^'^^.G^ ctgeac) leads to a non-convergent 
search for the minimal x^. Therefore, to pin down also the "Regge slope" parameters, we use a 
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Figure 10: Effective slope parameter W^, defined in Eq. (I92p . obtained from tlie S-PW model with dipole 
(1 -t/M2)-2 [a' = O.lS/GeV^) (dot-dashed) and exponential exp(t/26) (a' = 0) (dotted) t-dependence. 
Parameters where obtained by a NLO (MS) fit and are specified in Table HI Solid line is the fit (j90p of 
the HI Collaboration. 

three-step fitting procedure (fitting to DIS data being the first step). We first release the gluonic 
residue t-slope parameter, which yields for the exponential t-dependence 

E^*^'^ = 2.85 GeV-^ s'^'^^ = 0.71 , 5^ = 2.71 GeV-^ = 0.36 . (97) 

Now releasing the a'^^^ = = parameters (thick dashed hne in Fig. [TTl) . we find that they 
essentially stay at zero. Compared with our three-parameter DVCS fit ansatz (18^ (thin dashed), 
where fixed parameters were taken from J/"^ photoproduction, the gluonic t-slope only slightly 
increases (17%). Hence, we consider the two-gluon GPD models, obtained from DVCS and J/\E' 
photoproduction, as compatible. Moreover, the small values of the skewness parameters indicate 
that the exponential ansatz in the MS scheme is mostly compatible with the 1-PW model, as 
pointed out in Sect. 14.31 

For the dipole ansatz the parameter set after the second fitting step reads 

^^sea^)2 _ Q 54 Qgy2 ^ ^sea _ 5 24 , {M^f = 0.47 GcV^ , = 8.71 (98) 

with Ogpa^ = ttQ = 0.15/GeV^. We realize that the quark skewness and cut-off mass parameters 
only slightly change, cf. Table HI The gluon cut-off mass moderately decreases and the skewness 
parameter increases, compensating each other in the CFF and so the cross section is well described. 
Now releasing the a' parameters we essentially find that for the dipole ansatz (thick solid) they 
slightly reduce for quarks and more significantly for gluons: 

a^^^ = 0.15 ^ 0.12 , a'c = 0.15 ^ 0.08 . (99) 
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The quality of fits, compared to those shown in Table [3] for fixed gluon slope parameters, does not 
change. 
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Figure 11: Quark and gluon transverse distribution widths y (^^) with respect to for x = 10 ^ (a, 
b) and with respect to x for = lOGeV^ (c, d), obtained from NLO MS DVCS fits of S-PW model 
with dipole (sohd) and exponential (dashed) t-dependence. Thin lines correspond to parameter choice 
from (j83p or (j84p . while for thick lines first the gluon residual t-slope parameter or B'^ and then 
both Ogpa^ and parameters were released. 

Fig. [11] summarizes our NLO (MS) findings for the quark and gluon transverse width \l (6^), 
defined in Eq. flM|) . We note that for fixed gluon slope parameters (thin lines) the gluonic trans- 
verse width, shown in the right panel, would fully coincide with our LO curves. The only moderate 
difference appears in the dipole ansatz (thin solid), where the quark transverse width decreases 
from ^ 0.9 fm, quoted in Eqs. fl95|96l) . to ~ 0.8 fm. If we take the gluonic parameters from the 
three-step fit, we see that the gluonic transverse width for both the dipole (solid) and exponential 
(dashed) ansatz increases from 0.7 fm and 0.6 fm to ^ 0.85 fm and f« 0.65 fm, respectively. 
However, this moderate difference for the gluons affects only slightly the quark transverse width. 
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see the left panel. At NLO we observe already at the input scale that quark and gluon transverse 
width are mostly the same. The values we can quote 

r dipole 



0.75 - 0.80 



^ < > fm for ^ 

sea V ' 'G 10.60 — 0.651 [exponential 



ansatz . 



(100) 



are stable under evolution. Our results are fully robust with respect to scheme conventions or 
radiative NNLO corrections. 
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Figure 12: Quark (a) and gluon (b) transverse profile function (llOip for = 4GeV^ and x = 10 ^ is 
obtained from a six parameter DVCS fit, as in Fig. [TTl 



We conclude that within our GPD models, having a flexible skewness ratio and the same func- 
tional form of t-dependence for sea quarks and gluons, the resolution of the transverse distribution 
of partons is robust and that at small-x (sea) quarks and gluons have the same transverse width. 
For the exponential ansatz the DVCS result ^ 0.65 fm is compatible with the one quoted in 



Ref. 



271 ]. namely 0.65 fm. For the dipole ansatz, we find a larger value 



reflects the functional form we have employed. It is illustrated in Fig. 
the (normalized) transverse profile function 



0.8 fm, which simply 
that in impact space. 



(101) 



for dipole and exponential t-dependence mainly differ for distances larger than the disc radius of 
the proton, i.e., for b > 0.6 fm. Hence, the larger value of the transverse width for the dipole 
ansatz arises from the long range tail of the profile function, see the solid curve. Note that the 
model uncertainty in the extrapolation of the GPD to t = corresponds to the uncertainty in the 
long range tail. 



51 




Figure 13: cos(j) harmonic (dash-dot-dotted) of beam charge asymmetry (BCA) for E-PW model 
fit at LO as a function of azimuthal angle (p (a) ot of anomalous gravitomagnetic moment of sea 



quarks B^^^ = —B'^ (b), compared to preliminary HI data 10|]. Gray band is 1-a region of HI fit 



to pi COS0, pi = 0.17 ± 0.3 ± 0.5, with errors added in quadrature 

4.5 Is the anomalous gravitomagnetic moment accessible? 

Let us first adopt the "classical" Regge point of view in which the chromomagnetic "pomeron" is 
absent, i.e., we can neglect as above the OFF S. Hence, the BCA ([6]) is already predicted by the 
outcome of our cross section fits from the previous sections. A preliminary BCA measurement of 



the HI Collaboration with uncorrected acceptance effects has been reported in Ref. 10|| for the 
kinematics 

0.05GeV2 < |t| < lGeV^ {W) = 82GeV, and {Q'^) = 8GeV\ (102) 

To evaluate the BCA ([B]), we integrate over the t-interval and take the given mean values. In Fig.fT^ 
(a) we confront our LO prediction (solid curve) for the E-PW model (155]) with the preliminary 
BCA measurement. Having in mind that we neglected the cos(30) harmonics, induced by gluon 
transversity, and the presumable small twist-three cos(2</)) harmonics, we can safely state that our 
prediction is fully compatible with the preliminary HI measurement. 

As explained in Sect. [21 accessing the twist-two CFFs through the cos(</)) harmonic in the 
Fourier decomposition (ITO!) is in the considered kinematics relatively clean. We recall that this 
method can only diminish the contamination from other CFFs. The preliminary HI fit, 

pi = 0.17 ±0.3 ±0.5, (103) 

is shown in Fig. [T^ a) as error band and agrees well with our prediction (dashed curve), evaluated 
from the expression (j6]). The deviation of the dashed curve from the solid one is of kinematical 
origin and arises from the BH terms. In the absence of large genuine twist-three contributions and 
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gluon transversity or its twist-four contamination, all other harmonics in the Fourier decomposition 
(HU]) of the BCA are predicted to be small [56|: 

0<J>o<10-^ p2~-5■10^^ 0<ps<10-\ (104) 

where higher harmonics for i > 4 can be considered as negligible. We emphasize that a significant 
deviation from these numbers would indicate that genuine twist-three and gluon transversity 
(or twist-four corrections) are addressable in the small-xe kinematics, where it is assumed that 
electromagnetic radiative corrections will not alter the angular dependence of the BCA ([5]). 

To have a closer look to the parameter dependence of our prediction, we write the approxima- 
tion f lTT]) for t-integrated data (valid for ^ — tmin ~ ^B^p) 



Pi ~ n- 



with U = -0.05 GeV^ , to = -1 GeV^ 



W=S2 GcV 
Q2=8GcV2 



where n and h are two kinematical factors. Apart from the skewness effect, the size of pi is 
essentially governed by the "pomeron trajectory" 

p, oc cot(vra((t), {Q'))/2) ^ {a{{Q')) - l) , 

which for the mean values {-t) ^ 0.2 GeV^ and (Q^) = 8 GeV^ can be replaced by the intercept 
tt((Q^))- The justification for neglecting the a' parameter can also be read off from the last column 
in Table [21 which shows that both ansaetze, with a' = 0.15/GeV^ and a' = 0, yield almost the 
same value for the BCA. Hence, in our Regge-inspired GPD framework the BCA prediction arises 
to a great extent from the LO DIS fit. This is also accompanied by a large uncertainty, namely, 
a small absolute error of the DIS intercept a — 1 will induce a large relative one for the BCA. 
Certainly, it would not be appropriate to quote here the width of minimization curve as an 
error estimate, since our model simplifications, e.g., neglecting "Reggeon" contributions, might 
influence the fit result for a — 1, too. At NLO (NNLO) we observe that the BCA decreases to 
Pi = 0.9(0.12), see last column in Table H] ([6]), which corresponds to a 1-a deviation from the 
preliminary HI mean value fllOSp . The relative drop of the BCA size at NLO, compared to LO, 
seems indeed to be correlated with that of the intercept, cf. corresponding columns in Tables [2] 
andU 

We mentioned in Sect. [2] that a measurement of 



Pi ~ n- 



W=82 GoV 
Q2=8GcV2 



h+ei::dt{\w-jij,\s\ 

supplemented by the DVCS cross section measurement, should in principle allow to separate the 
CFFs "H and £. It is obvious from the formula that a negative real part in S will more strongly 
influence our BCA "prediction" than a positive one. 
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To quantify the sensitivity of the BCA in dependence on B^'^'^, we adopt the GPD E model 
fl58p and rely on the scenario B'^ = —B^'^'^. We vary B^'^^ at the input scale from —0.4 ■ ■ ■ 0.4, which 
corresponds to a variation of 

-0.13 < J^'^^ < 0.28 (0.43 >J^> 0.02) with = 0.2. 

In the right panel of Fig. [13] we show the dependence of the BCA (at = 0) on the B^^^ parameter 
for the S-PW model at LO accuracy. We refitted for each given value of ^Bsea the model parameters 
of the GPD H so that the DVCS cross section is described. As one realizes, pi depends only slightly 
on the parameter B^^'^; its relative change is 

~ -15%--- + 5% for -B'"^ = -0.4 ■ ■ ■ 0.4. 

We obtain rather similar findings for a B*^ = scenario and within the nl-PW model. Contrarily 
to our hopes, we conclude that both the model uncertainties of "H and the experimental errors are 
too large to obtain a bound for B^'^'^. 

4.6 Lessons from fits 

We would like now to summarize the lessons from our fits, presented in the previous sections. We 
recall that our DVCS description relied on Regge-inspired GPD models, set up at a (low) input 
scale, and the coUinear factorization framework at leading twist-two in LO and beyond. Both the 
factorization and renormalization scale were set equal to the virtuality of the incoming photon. 
Concerning strategies for fitting to the unpolarized DIS structure function and DVCS cross section 
measurements we have seen that 

• GPD models with a flexible skewness ratio describe well the small-xe DVCS cross section 
measurements at LO and beyond, 

• fitting parameters obtained from simultaneous (DVCS+DIS) and two-step (first DIS, than 
DVCS) fits are rather similar, 

• and that the partonic shrinkage effect can be addressed by a three stage fitting strategy. 

Our Regge-inspired modelling within the coUinear factorization framework, including the corre- 
sponding approximations, is consistent with experimental measurements: 

• small "Regge slope" parameters a' < 0.25/GeV^ are favored and large ones for "pomeron" 
related quarks, e.g., a' ~ 1/GeV^, are excluded. 
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• the scale dependence of both the residual t-dependence and "Regge slope" are compatible 
with pQCD, 

• and the real part of the amplitude, following from DIS and DVCS cross section fits, predicts 
the preliminary BCA measurement. 

Relying on the assumption that the CFF "H is dominant and that higher twist contributions 
are unimportant, we conclude from our fits in Sect. 14.21 that realistic GPD models at LO have the 
following properties: 

• quark GPD models possess a skewness ratio r^*^^ ^ 1 

• and to ensure this for a large lever arm, gluon GPDs have a skewness ratio r'^ < 1. 

We emphasize that our finding r'^'^^ ^ 1 was predicted by the aligned-jet model and it is realized 
in an RDDA for W^'^ ^ 1. However, all popular gluonic GPD models, possessing f« 1, are 
incompatible with the small-xe DVCS data (and presumably also in a pure LO description of 
data on hard vector meson electroproduction, see, e.g., the dashed curves in Fig. 5 of Ref. j46| 
that refer to the coUinear factorization prediction and are included there for illustration). The 
negative skewness ratio of our gluonic GPD models at the input scale and the feature that it 
approaches the conformal ratio for growing indicate that the parameterization of the skewness 
function is rigid. We can safely state that 
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961] that the skewness ratio is equal to its conformal value is 



• the small-x claim of Refs. 
ruled out at LO 

• and a flexible parameterization of evolution has not been achieved so far. 

The reader might be surprised by this second conclusion, which contradicts opinions and model- 
dependent findings in the literature. Nevertheless, we are not aware of a distinct investigation 
that addresses the evolution of the residue function in dependence of the initial condition, i.e., the 
skewness function in momentum fraction space, the series of conformal PWs in J space, the series 
of forward-like functions Qp{z) in the "dual model", or the series of "conformal sibling poles" in 
Mellin space. 

Including radiative corrections, we observed in Sect. 14.31 that NLO GPDs qualitatively differ 
from the LO ones: 

• the quark skewness ratio in the MS scheme can be larger than the conformal one, 

• the quark skewness ratio in the CS scheme matches the conformal one. 
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• the scheme dependence at NLO is absorbed by a sizeable GPD reparameterization, 

• and the gluonic skewness ratio approaches in both schemes the conformal one. 

The large reparameterization effects, in particular for quarks, might be a consequence of the 
remaining rigidity of the skewness function in our flexible models. We expect that this can be 
overcome by taking into account two effective non-leading S0(3)-partial waves rather than one. 

We would also like to emphasize that the scheme dependence at NLO, we pursued in our 
studies, is entirely related to a 'non-diagonal' rotation, which shows that 



the small-x claim 
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96| of conformal skewness ratio can be sustained beyond LO only in 



special schemes, however, presumably not in the MS one. 

We add that our MS fit at NLO for GPD models with an exponential t-dependence is consistent 
with the LO conformal skewness ratio. In our opinion this fact does not necessarily support the 



"logic" of the small-x claim, since it relies on tree-level conformal symmetry 9^, |96| that is beyond 
LO not explicitly implemented in the MS scheme. Hence, the "predicted" skewness ratio in the 
MS scheme at NLO differs from the LO one. 

The stability of the perturbative approach can be presently studied only in the CS scheme, 
and there 

• the inclusion of NNLO corrections yields only a small change of fitting parameters. 

Thereby, the skewness ratio in the CS scheme is essentially given by the conformal ratio ( 128|1 . 
Thus, in this scheme beyond LO one might be tempted to conjecture that a GPD is indeed tied to 
the corresponding PDF by conformal mapping. One might wonder whether it is accidental that 



this hypothetical 'holographical principle' is realized within the group S0(2,l) [64 1. 

The functional form of t-dependence cannot be pinned down from present DVCS data. How- 
ever, we can definitely state that: 

• "Regge slope" parameters are small at the input scale, however, do not necessarily vanish, 

• non-vanishing "Regge slope" parameters rapidly decrease with growing scale, 

• both residual dipole and exponential t-dependence are compatible with present DVCS data, 

• within our models, there is a cross-talk between skewness and t-dependence at LO, 

• and beyond LO the skewness and t-dependence start to decouple in our flexible models. 
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In a three-step fitting strategy we also observed a correlation between the functional form of the 
residue function and the partonic shrinkage effect. Namely, an exponential t-dependence of the 
former is quite consistent with a zero shrinkage effect, i.e., a'^^.^ = o'q = 0. 

Our partonic interpretation of the transverse degrees of freedom, given in Sect. 14. 4[ is rather 
robust. Relying on the decoupling of t- and skewness-dependence, we have found for dipole 
(exponential) ansatz that 

• sea quarks and gluons have roughly the same transverse width ~ 0.75(0.65) fm, 

• the transverse width is rather stable under evolution, 

• and a possible partonic shrinkage effect at a lower scale is rapidly washed out with growing 
scale. 

Unfortunately, we saw in Sect. 14.51 that both theoretical uncertainties and experimental errors 
do not allow to address 

• the chromomagnetic "pomeron" or angular momentum of sea quarks, 

• and gluon transversity. 



5 Small- ccb fit results as input for "dispersion relation" 
fits 

The reader might wonder why have we repeatedly stressed that the claim [qgJ that at small-x 
skewness ratio has a conformal value (128|) is excluded at LO, when this value seems to be realized 
in the CS scheme beyond LO. However, this is entirely related to the ambiguities present in our 
factorization conventions and partonic interpretations. For instance, in a DVCS scheme, where 
the gluons are not resolved in the hard-scattering subprocess, above LO findings at the input scale 
would not be essentially altered by radiative corrections, which only modify the pQCD evolution 
predictions. 

Such an interpretation is favored among model builders and also in the GPD phenomenology 
of fixed target kinematics. We now utilize our LO analysis of the small-xe data for a "dispersion 
relation" fit in fixed target kinematics, where real part of OFF is taken to be determined by 
imaginary part and subtraction constant, as described in Sect. [2] (see also Ref. 6J] for more 
sophisticated strategies and Ref. |l39| for an alternative approach). This substantially reduces 
the model uncertainties in common GPD model description. For fixed target kinematics, where 
the lever arm is rather limited, one may additionally rely on the so-called scaling hypothesis. 
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i.e., on the assumption that the GPD does not evolve under the change of the photon virtuahty. 
The primary goal of such fits is to reveal the shape of the dominant GPD H on its cross-over line 
(?7 = x) from DVCS measurements on unpolarized proton target. 

The framework, as described in Sect. [2] for small x-q can be easily adapted. It is beyond 
the scope of this paper to present a detailed description of our fits; however, we would like to 
demonstrate that GPD phenomenology for DVCS can be straightforwardly set up. We consider 
now three active quarks and write the partonic decomposition of '^mH as 



e, t) + e, t) + Ih^'^^ t) 



(105) 



taken at an "input scale" of 
DD representation ( l23l) . 

F{x, X, t) = 



2 GeV^. We model the GPD on the cross-over line using the 



^ + X Jq \1 + X 1 + X 

and take the t-dependence from the quark spectator model 
functional form: 



2u 



X 



H{x, X, t) 



n r 



2x 



1 + X \1 + X 



1 — X 



X 



1 



V l+xhP) 



(106) 

This suggests the following 

(107) 



Here n is the residue normalization of the PDF, which can be taken from PDF fits, r is the 
skewness ratio at small x, a{t) is the "Regge trajectory", which can be borrowed from Regge 
phenomenology, h controls the large x behavior, which according to counting rules 140| should 
be different from that of PDFs, and both p and M control the t-dependence. The functional 
dependence on t and x in the t-dependent part of the ansatz is specifically motivated by the 
spectator model 88|. The small-x behavior of the sea quarks is taken from our small-xe DVCS fit 
at LO, evolved backwards. The corresponding GPD is described by the model (11071) within the 
parameters 



d'^'^t) = 1.13 + 0.15 t/GeV^ 



n 



1.5 



(M' 



/rsca \ 2 



0.5 GeV' 



2. (m 



For valence quarks, in accordance with our Regge-inspired modelling, whi ch im plies that we 
take the p and u trajectories, and standard PDF parameterizations, e.g., Ref. 1301 ]. we fix quark 
parameters to be: 



a™Xt) = 0.43 + 0.85 t/GeV' 



n 



, val 



1.35 , (M 



vah2 



0.64 GeV^ 



^val 



P 



(109) 



i.e., we take a monopole ansatz with fixed cut-off mass for the residual t-dependence. Let us argue 
that our model for the GPD on the cross-over line is generically consistent with some full valence- 
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quark GPD model that satisfies the form factor sum rulqlj- Our choice p™' = 1 is motivated by 
the dimensional counting rule for the Dirac form factor, which states that this form factor behaves 
as (— for — t — 7- oo. Loosely spoken, the missing power here would arise in our GPD model 
from the x- integration for fixed 77, e.g., rj = 0, providing us another monopole with cut-off mass 



a 



val 



,/val 



M 



val 



0.8 GeV. 



This choice is consistent with the characteristic nucleon scale and our ansatz is roughly asso- 
ciated with an effective dipole parameterization for the Dirac form factor with a cut-off mass 
of M™' ^ 0.8 GeV. Such an effective parameterization of the Dirac form factor is for small 
— t < 0.5 GeV^ region, on which we are interested, at the 20% level consistent with the stan- 
dard dipole parameterization, used above in Eq. 



Furthermore, we take into account the pion-pole contribution 149| . displayed in Eq. (153) 
of Ref. {551, which, however, is a non-dominant contribution in DVCS. More important is the 
subtraction constant f lT^ 

C 



Cit) 



2 ' 



;iio) 



where the normalization C and dipole cut-off mass are fitting parameters. According to the 
"dispersion relation" f|T2|) . we also include this constant in the non-dominant OFF £ = C{t), 
however, set its imaginary part as well as the CFF "H to zero. 
The five fitting parameters, 



„val 



ival 



c, M^"^ 



control the large- and small-x behavior of the imaginary and the real part of T-L. Note that by 
a "dispersion relation" fitting procedure one accesses even the large, experimentally inaccessible 
X region. We also recall that the subtraction constant (1121) can be eval uated from the so-called 
D-teim and it should not be confused with a so-called J = (fixed) pole 
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1421 in the context 



of merging parton model and Regge phenomenology, see footnote [31 the discussion below Eq. 



in Sect. 13.21 and for more details Refs. 62l. |64|. 



In the fitting procedure we utilize the formula set of Ref. 56|] , which are based on 1 / 
approximation for the squared amplitude. Unfortunately, this approximation is inappropriate for 
JLAB kinematics; still, it can be improved in the unpolarized target case by a 'hot fix', taken from 



^^Note that in popular Reggc-inspired GPD models, used in DVCS phenomenology, the form factor constraint, 
which is already lost on generic level, is not implemented. Unfortunately, such models are sometimes employed to 
"constrain" quark angular momentum. 
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Figure 14: Description of data by global model fits to HERMES |145l | and CLAS [146l | data (squares, 
dashed line) and additionally including the HALL A measurements (circles, solid line): (a) BCA A^^^ 

from HERMES versus t (diamonds), (b) BSA A{(f> = 7t/2) from CLAS versus for t 0.15 GeV^ 

(upside triangles) and t ~ — 0.3GeV'^ (triangles), (c) helicity-dependent and (d) unpolarized differential 
cross section from HALL A for = 0.36, —t = 0.33 GeV^, and = 2.3 GeV^ (diamonds). Note that 
in (d) the cos(2(/>) harmonic of the interference term is not included. 



are the BCA from HERMES 



144 



a refined evaluation for a scala r target 143l | . The available fixed target data for unpolariz ed p roton 



145| . the beam spin asymmet ry (B SA) from CLAS 1461 ]. and 

' 147I . The covered kinematics in 



the electroproduction cross section measurements from HALL A 
^(= x) is depicted below in panel (a) of Fig. [T5J 

Within the hypothesis that "H is the dominant CFF we can describe HERMES and CLAS 
data, where in the later case only small —t <^ data were taken into account. Altogether we 
included 36 data points and found with 

xVd.o.f. ^ 32/31 

an acceptable fit, see squares in panels (a) and (b) of Fig . [TH The fit also predicts a preliminary 



BSA measurement from the HERMES Collaboration 



148l ] (not shown) as well as, to some extent. 
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the helicity-dependent cross sectio measurement from the Hall A Collaboration, however, not 
the unpolarized (helicity-independent) one at —t = 0.33 GeV^, displayed by dashed curves in 
panels (c) and (d) of Fig. [HI respectively. It provides the following parameter set for our GPD 
model: 



3.09, 



,val 



0.95, b 



val 



0.45 



C = 0.24, M' 



sub 



0.5 GeV. 



The b val ues are smaller than the corresponding /3 ones for PDFs, which is in accordance with 



Ref. 



1401. The skewness ratio is smaller than the conformal one, which is in this case ~ 1.2. 



However, one should bear in mind that we did not include "Reggeon exchanges" in the flavor 



smg 



150 



et se ctor. The para meter of C is here smaller than suggested by chiral quark soliton model 
calculations and the t-dependence of the subtraction constant is rather 



153j or lattice 
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steep. 



The hypothesis of the ff-domi nance is within our model ( I107p no longer valid for the cross 



147| . performed at large xb = 0.36, = 2.3 GeV^, and for 



section measurements of HALL A 
four t values, 

-t = {0.17 GeV^ 0.23 GeV^ 0.28 GeV^ 0.33 GeV^}. 

To describe the full set of these HALL A data, it is required that the real part of the DVCS 
amplitude steeply varies in this t interval and becomes larg e at -t = 0.33 GeVl Both of these 
features are not implemented in common GPD models 154l Il55| . We will now explore whether 
these features of the real part arise from the inclusion of the GPD H. To be specific, we take a 
model for the CFF 



IT 



4 1 
2- + - 
9 9 



n 



2e 



-a™'(t) 



1 



3/2 



1+5 (M^'>-i)2 

where for simplicity we employ the parameter set (11091) and fixed b = 3/2. Hence, only the 
normalization n appears as a new fitting parameter. If we take this parameter of order one, our 
ansatz is to some extent consistent with standard GPD models. 

In our second global fit we take the DVCS data set from above and include four values of the 
azimuthal angle asymmetry 



aw cos[(p)- 



dM<l^, ...) 

aw -—— 

dcpdx-QdtdQ^ 



iQ d(f)dxBdtdQ'^ / jQ 

which are obtained from the four unpolarized cross section measurements of the HALL A Col- 
laboration. Here, the integral measure dw, defined in Eq. (103) of Ref. 5^, compensates the 



^° We adopt here the terminology of HALL A [147| , where this observable is defined as half of the cross section 
difference for positive and negative electron beam helicity. 
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dependence of the Bethe-Heitler propagator product. In other words, in our fit only the relative 
changes of the unpolarized cross section with respect to a cos modulation, dominated on ampli- 
tude level by twist-two CFFs, is essential. Within such a strategy the absolute normalization of 
the cross section is irrelevant. 

Compared to the previous global DVCS fit, this provides us four data points and one degree 
of freedom (normalization of "H) more. The quality of fit 

xVd.o.f. ^ 33/34 

is good and the parameters for the GPD H are 

b'^^ = A.6, r™' = l.ll, 6™i = 2.4, C = 6.0 , M™^ = 1.5GeV. (112) 

The value of the subtraction constant agrees now qualitatively with the original chiral quark 



soliton model estimate |l50l |. quoted as C ~ 5. For the normalization of "H we find the value 
n ^ 3. This corresponds to a huge skewness effect ~ 5. We consider this finding as an effective 
parameterization of some GPD contribution that reveals our lack of partonic understanding of the 
unpolarized cross section measurement of the HALL A Collaboration. We also like to stress that 
the real part of the DVCS amplitude at larger values of xb is an intricate quantity that cannot 
be simply related to a specific GPD property. Note also that a closer look reveals that essentially 
only the HALL A data point at —t = 0.33 GeV^ forces us to explore the scenario of non-standard 
GPD modelling. 

In this w ay, o ur model (11121) describ es p resent DVCS data, including the helicity-dependent 



cross section |l47l | and preliminary BSA |148l | data which have not been employed in the fits. This 
is exemplified by the circles and the solid curves in Fig. [TH Note that in panel (d) the deviation of 
our fit result from the data points might be attributed to a cos(20) harmonic in the interference 
term. However, the extraction of GPDs from a DVCS measurement is highly nontrivial and we 
relied here on model assumptions and employed a simple least square fit. Thus, we like to spell 
out a clear warning, namely, the result we presented does not exclude a successful description of 
the data within a rather different set of parameters: 

• a literal interpretation of these (first) model-dependent fit results is not appropriate. 

Nevertheless, the parameters (11121) look rather reasonable from the generic point of view, as stated 
above, and they might be used to set up an H GPD model in any favored representation. 

The resulting GPDs on the cross-over line are shown in the left panel of Fig. flT5]) . It is 
obvious that the two different fits provide results which are rather different in the large-x region, 
but approach each other at smaller x. Essentially, large difference for JLAB kinematics simply 
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Figure 15: (a) Global GPD H fits (fTTT]l and (fTT2]) at t = -0.3 GeV^ (thick) and t = (thin) are 
displayed as dashed and solid lines, respectively, (b) Prediction of the BCSA asymmetry ()113p for 
COMPASS kinematics (S^ = 160 GeV, = 2GeY^, t = -0.2 GeV^) versus v? = vr - 

indicates that BSA measurements can be described with two qualitatively different GPD model 
scenarios, where the dashed curve represents the common GPD models with a, let us say, moderate 
DVCS amplitude, while the solid curve belongs to GPD models that provide an enhanced DVCS 
amplitude, here induced by H. Whether this is a realistic scenario is a open problem, which should 
be addressed in future studies. 

We can now employ our model fits to deliver a prediction for the COMPASS experiment. In 
such a fixed target experiment one would scatter positively or negatively charged muons with 
helicity —1/2 and +1/2 on a proton target. The preferred observable is the beam charge-spin 
asymmetry (BCSA), 



A 



BCSA 



(0) 



;ii3) 



which is essentially related to the real part of the interference term. It is therefore rather sensitive 
to details of the spectral function Qm'H. We recall that the sign of the BCA and, thus, also of 
the real part of the CFF "H changes somewhere between HERMES and HI /ZEUS kinematics. 
Knowing the position of this zero in dependence on the kinematical variables would be crucial 
for pinning down GPDs. For the mean values at COMPASS, we expect a sizeable asymmetry, 
the sign of which is fixed by valence-like quarks (or "Reggeon exchange"), where the subtraction 
constant (fixed pole) plays a role, too. We display our predictions resulting from our two fits 
for the BCSA f lllSp in the right panel of Fig. [15] versus the azimuthal angle ip = n — (p (i.e., 
within a so-called Trento convention). Therefore, the COMPASS kinematics x ~ [10^^, 10^^] is 
certainly well suited to explore, in the LO DVCS interpretation, the transition area between sea 
and valence-like quarks domination. In the NLO interpretation this translates into the hope of 
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exploring the interplay between quarks and gluons. 

6 Summary and conclusions 

In this paper we demonstrated that for the description of small-xe DVCS data the choice of 
representation, used to set up GPD models, does not matter too much. We employed different 
"languages" and found that, although explicit transformation formulae are only partly available, 
one can easily translate main results. The problem remains always the same, namely, to find a 
flexible parameterization of the skewness function and the skewness ratio. As a side remark, we 
explained that the small- a: claim, stating that the GPD is rigidly tied to the corresponding PDF 
by the conformal skewness ratio, is based on unjustified mathematical assumptions. Moreover, 
we demonstrated in detail that this claim is ruled out at LO. We also revealed that common 
gluon GPD models (either RDDA, t-decorated PDF, leading PW or minimalist "dual" model) 
possess the conformal ratio. We showed then that this is the very feature which is responsible 
for the failure of previous attempts to describe DVCS at LO, even if the quark GPD model is 
flexible. The "workaround recipe" to take NLO gluon PDFs in an LO description fails for our PDF 
parameterizations in DVCS. It also creates an obvious inconsistency because the same amplitude 
is described differently in DVCS than in its DIS limit. 

We introduced two flexible GPD parameterizations with respect to the skewness ratio, namely, 
one model was set up by adding a next-leading SO (3) PW and the other one by a model-dependent 
resummation of SO (3) PWs. Thereby, we relied on the simplest assumptions, e.g., (almost) 
decoupled skewness and t-dependence, and dressed the gluon GPD with the t-dependence from 
the J/\l/ production analysis. However, it turned out that these two models are rather similar and 
possess a rigid parametrization of the skewness function. Nevertheless, within these Regge-inspired 
models we could describe small-xe DVCS cross section data from the Hl/ZEUS Collaborations at 
LO and beyond. Thereby, the gluonic LO GPD on its cross-over line becomes negative at lower 
values of and, moreover, we had a rather large reparameterization of the skewness effect at 
NLO for both the quark and gluon GPDs. Both findings are not expected and we consider them 
as model artifacts. So we conclude that the effective nl-PW, or the minimal "dual" , or any similar 
model is not flexible enough to control the evolution over a large lever arm in Q^, and so also 
the control of the skewness ratio at larger values of is lost. Also, the fact that at LO the 
skewness and t-dependence are still correlated in our model fits reveals that an effective nl-S0(3) 
PW parameterization still suffers from rigidity. 

We have demonstrated that at NLO the abovementioned small-x claim within the HERA- 
II DVCS data also does not necessarily hold true. Within a dipole t-ansatz, however, we have 
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confirmed our previous results that claim is to a great extent valid in the CS scheme, at NLO 
and NNLO. Thereby, we observed very small NNLO corrections, namely, below 1%. This might 
indicate that the 'holographic' principle, which ties GPD and PDF, arises from a broken S0(2,l) 
symmetry. We stress that we do not consider such 'holographic' principle reliable enough to 
employ it in a fitting procedure, e.g., to reveal the t-dependence of GPDs. 

We utilized our flexible models to deliver the transverse distribution of quarks and gluons and 
found that they are quite robust with respect to the change of perturbative order and scheme 
conventions. We only observed a slight difference between LO and beyond LO results, which is 
attributed to the different fitting strategies and the remaining rigidity of our models. The results 
moderately differ from our previous findings, where we employed the 1-PW model (conforming 
to the small-x claim). In contrast to those findings, we observed now that within our flexible 
GPD models there is no significant difference between quark and gluon transverse width. The 
perturbative prediction that a partonic shrinkage effect should be washed out at larger scale 

~ 10 GeV^ is in agreement with experimental findings, however, its existence at the lower scale 
cannot be excluded. We employed a three-step fitting procedure to pin down the t-dependence of 
the GPD. Within an exponential t-dependence GPD model there is no shrinkage effect and the 
partonic transverse width is about ~ 0.6 — 0.65 fm. This value is compatible with the one extracted 
from the J/\E' photoproduction and with the radius of the proton disc. A dipole t-dependence 
ansatz is accompanied with a small shrinkage effect and the partonic transverse width is now 
0.75 — 0.8 fm. The quoted difference of widths arises essentially from the different extrapolation 
of the measured t-interval to t = and essentially provides an uncertainty in the long range tail 
of the profile function in impact space. 

Based on the DVCS cross section fits, we refined a previous model prediction for the beam 
charge asymmetry. Our prediction is compatible with preliminary measurements from the HI 
Collaboration, which supports our Regge-inspired GPD modelling. Unfortunately, both theoretical 
uncertainties and experimental errors do not allow us to access the chromomagnetic "pomeron" 
and so neither the anomalous gravitomagnetic moment (or angular momentum) of sea quarks. 
An immense reduction of both sources of error is needed before one might be able to address 
these questions. An improvement of the experimental result might lead to some clarification of 
the presence of a third azimuthal angular harmonic, induced by gluon transversity. 

Based on the double distribution representation, we have also built a simple model for the GPD 
at the cross-over line, where we implemented the small-x behavior, extracted from Hl/ZEUS data. 
Using this model and relying on the scaling hypothesis, we presented a first "dispersion relation" 
fit to observables for fixed target DVCS experiments on unpolarized proton target. Assuming 
some unexpected properties of if, we were able to describe all available small —t data in a first 
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global fit. The parameters we found match our generic expectations coming from Regge behavior, 
large-x counting, and quark soliton model estimates. Utilizing this GPD as an input, we predict 
beam charge-spin asymmetry as measurable at COMPASS. Thereby, it becomes obvious that 
COMPASS is needed to reveal the GPD H further. This experiment might also be important for 
pinning down the transverse distribution of partons in the transition region between "pomeron" 
dominance and "Reggeon" behavior. 

Let us finally emphasize that one can have different partonic interpretations of our findings. 
Namely, one can either state that gluons are very important in DVCS and should be perturba- 
tively resolved or one sticks to the quark picture (DVCS scheme) in which gluons only enter the 
evolution equation. This is a matter of taste that so far simply provides two qualitatively different 
GPD parameterizations and also points in two different directions for future DVCS studies. We 
emphasize that in any case an LO description is a convenient starting point for the phenomenolog- 
ical description of the growing amount of experimental data. It offers the possibility to pin down 
GPD models within rather straightforward fitting strategies. Such models are phenomenologically 
valuable for making contact with dynamical models and lattice simulations, since in the latter 
case matching of lattice and perturbative renormalization is for GPD moments only done at LO. 
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Note added: After our manu scrip t was finalized, we have noticed the new DVCS measurements 
from the ZEUS Collaboration 156| . We have convinced ourselves that these data will only mildly 



influence our fits and, thus, we do not update here our results. In other words, we might view our 
findings as predictions for the new ZEUS data. After the experimental DVCS analysis of HERA 
II run will be finished, we plan to include them in a forthcoming update. 
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